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^ ■ ABSTRACT 

We present multi-epoch, near-infrared imaging of the binary T5.0+T5.5 dwarf 

2MASS J15344984-2952274AB obtained with the Keck laser guide star adaptive optics 

system. Our Keck data achieve sub-milliarcsecond relative astrometry and combined 

1^ ' with an extensive (re-)analysis of archival HST imaging, the total dataset spans ~50% 

I ■ the orbital period. We use a Markov Chain Monte Carlo analysis to determine an orbital 

\^ . period of IS.lj^^'g yr and a semi-major axis of 2.3I'^q'2 ^U. We measure a total mass of 

c/3 , 0.056 lb 0.003 Mq (59 lb 3 Mjup), where the largest uncertainty arises from the parallax. 

This is the first field binary for which both components are confirmed to be substellar. 



(N 



X 



This is also the coolest and lowest mass binary with a dynamical mass determination to 



^ ' date. Using evolutionary models and accounting for the measurement covariances, we 

-^ ' derive an age of 0.78ib0.09 Gyr and a mass ratio of 0.936^Q'QQg. The relatively youthful 

CN ' age is consistent with the low tangential velocity of this system. For the individual com- 

P : ponents, we find T^s = 1028 ib 17 K and 978 ±17 K, log(g) = 4.91 ±0.04 and 4.87 ±0.04 

^ : (cgs), and masses of 0.0287 ± 0.0016 Mq (30.1 ± 1.7 Mj^p) and 0.0269 ± 0.0016 M© 

OO ! (28. 2 ± 1.7 Mjup). These precise values generally agree with previous studies of T dwarfs 
O 
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and affirm current theoretical models. However, (1) the temperatures are about 100 K 
cooler than derived for similar field objects and suggest that the representative ages of 
field brown dwarfs may be overestimated. Similarly, (2) the H-R diagram positions are 
discrepant with current model predictions and taken at face value would overestimate 
the masses. While this may arise from large errors in the luminosities and/or radii 
predicted by evolutionary models, the likely cause is a modest («100 K) overestimate 
in temperature of T dwarfs determined from model atmospheres. We elucidate future 
tests of theory as the sample of substellar dynamical masses increases. In particular, we 
suggest that field brown dwarf binaries with dynamical masses ("mass benchmarks") 
can serve as reference points for Tgg and \og{g) and thereby constrain ultracool atmo- 
sphere models, as good as or even better than single brown dwarfs with age estimates 
("age benchmarks"). 



Subject headings: binaries: general, close 
techniques: high angular resolution 



stars: brown dwarfs — infrared: stars 



1. Introduction 

Over about the past decade, the parameter space of traditional stellar astrophysics has been 
greatly expanded with the discovery and characterization of brown dwarfs, objects that for most 
of their lifetimes are colder and less luminous than main-sequence stars. Despite ample progress in 
finding and characterizing brown dwarfs, very few direct measurements of their physical properties 
have been made so far. In particular, dynamical masses for brown dwarfs are sorely needed to test 
the theoretical models over a wide range of parameter space. In comparison to the >100 binary 
stars with direct mass determinations, dynamical masses have been measured for only a handful of 
objects clearly below the stellar/substellar boundary: 



1. the secondary component of the M8.5+M9 visual bin ary Gl 569Bab, which itself is a com- 



panion to a young (r:! 1 00-300 Myr) field M2 5 dw arf (JMartin et al.l l200d : iLane et al. 



Kenworthv et al. 



2001 



Zapatero Osorio et al.ll2004l ): 



2001 



2. +3. the very young (^My r) eclipsing M6.5-I-M6.5 binary brown dwarf 2MASS J05352184-0546085 
(jStassun et al.l 120061 ): 



4. the secondary component of the GJ 802AB system, where the secondar y has an estimated 
spectral type of L5 and the priinary is a M5.5+M5.5 spectroscopic binary ([Pravdo et al.ll2006l : 



Llovd et al.l 12009 : llreland et al.ll2008l ^ 



In addition, GJ 569Ba itself may be a unresolved binary brown dwarf (jSimon et al.ll2006l ). and the 
secondary component of the L0+L1.5 binary 2MASSW J0746425-h2000321AB, which appears to 



be an old (>1 Gyr) field system, has a mass near the stellar/substellar boundary (JReid et al.ll200ll : 
Bouv et al.l[20oi : lGizis k Reidlbood l. 

About 100 ultracool visual binaries are known£| found with high angula r resolution imag- 
ing surveys conducted by Hubble Space Telescop e {HST; e.g.. iReid et al.ll200ll : iBouy et alJ 120031 : 



Gizis et al. 



2OO3I: iBurgasser et al.l l2003bl. l2006d^ and ground-based adaptive optic s (AO) imag- 



Close et al. 



200c 



31: IBurgasser et al.ll2005al : iLiu et al.ll2006bl : iReid et al.ll2008l : Liu et al., in 



ing (e.g 

prep). Only ~10% of these are binaries composed of the coldest class of brown dwarf, the T dwarfs. 
T dwarfs are distinguished by their very red optical colors arising from pressure-broadened alkali 
resonance lines and very blue near-infrared (IR) colors from stroiig CH4, H2O, and c ollision-induced 
H2 absorption (e.g. lOppenheimer et al.lll995l : lGeballe et al.ll2002l : lKirkpatricldl2005l ). These are the 
lowest luminosity and coolest objects directly detected outside of our solar system, with bolometric 
of < 1O~^'^L0 and effective temperat u res (Tpff) of 



luminosities (Lh, 



2004 



qU 



Golimowski et al. 



Delorme et al. 



2004a 



700-1300 Kfe.g. 



Burgasser et al.l l2006al : iLeggett et al.l l2007al : I Warren et al. 



Vrba et a 



2007 



20081 ). As such, analyzing their physical properties is an important pathway to 



understanding the properties of gas-giant extrasolar planets. 

The subject of this paper is the T dwarf 2MASS J1534 4984-2952274AB, her e inafter 2MASS J1534-2952AB, 
which has an integrated-light infrared spectral type of T5 ( Burgasser et al.ll2002l . l2006b)o This ob- 



ject w as first resolved as a 0.065" binary in August 2000 in HST/WFPC2 imaging ([Burgasser et al. 



2003bl ) . Among known visual ultracool binaries i n the field, this system has the shortest estimated 
orbital period, only 4 years (see compilation in IBurgasser et al.ll2006dl^. In co mbination with its 
very high quality parallax measurement of 73.6 ± 1.2 mas (JTinney et al.ll2003l ) and the fact that 
its two components are nearly equal-magnitude (indicating nearly equal masses), this system is a 
very appealing object for astrometric monitoring. 

Laser guide star (LGS) AO provides a powerful tool for high angular resolution studies of brown 
dwarf binaries. Through resonant scattering off the sodium layer at ~90 km altitude in the Earth's 
atmosphere, sodium LGS s ystems create an artificial star bright enough to serve as a wayefron t 
reference for AO correction (JFoy &: Labeyridll985l : iThompson Sz Gardnerlll987l : iHapper et al.lll994l ) . 
Thus, most of the sky can be made accessible to near diffraction-limited IR imaging from the 
largest existing ground-based telescopes. We have previously used Keck LGS AO to discover that 
the nearby L dwarf Kelu-1 is a binary system (iLiu fc Leggettll2005l ) and to identify the novel L-l-T 
binary SDSS J1534+1615AB (|Liu et al.l l2006bl ) . In regards to dynamical mass determinations, 
ground-based telescopes equipped with LGS AO can provide the necessary long-term platforms 
for synoptic monitoring of visual binaries, especially where the required amount of observing time 



^We follow the popular convention that "ultracool" refers to objects of (integrated-light) spectral type M6 or later. 
Burgasser et all (|2006bl ) report integrated-light spectral types of both T5 and T5.5 for 2MASS J1534-2952AB 



based on a spectrum obtained with the CTIO/OSIRIS instrument. Examination of their spectrum by us finds that 
T5 is the correct typing. Also, a new near-IR spectrum obtained with the IRTF/SpeX spectrograph confirms an 
integrated-light spectral type of T5 (ij |2.3p . 



at each epoch is relatively modest but many epochs are needed, in contrast to HST where target 
acquistion can be slow and monitoring a populous sample over many epochs is quite telescope 
time-intensive. 

We present here the results of multi-epoch imaging of 2MASS J1534— 2952AB, observed as 
part of our ongoing high angular resolution study of ultracool binaries using LGS AO. § 2 presents 
our Keck LGS AO observations and (re-)analysis of archival ifST imaging. § 3 presents the resolved 
photometric properties of the binary and fitting of the orbit using a Markov Chain Monte Carlo 
method. § 4 compares the resulting total mass against evolutionary models, and § 5 summarizes 
our findings. Those readers interested solely in the results can focus on § 4 and § 5. 



2. Observations 



2.1. Keck LGS AO 



We imaged 2MASS J 1534- 2952 AB from 2005- 2008 using the sodium LGS AO svstem of th e 
10-meter Keck II Telescope on Mauna Kea, Hawaii (IWizinowich et al.ll2006l : Ivan Dam et al.ll2006l ). 
Conditions were photometric for all the runs. We used the facility IR camera NIRC2 with its 
narrow field-of-view camera, which produces a 10.2" x 10.2" field of view. Setup times for the 
telescope to slew to the science targets and for the LGS AO system to be fully operational ranged 
from 7-20 min, with an average of 12 min (e.g. lLiijl2006l ). 



The LGS provided the wavefront reference source for AO correction, with the exception of 
tip-tilt motion. The LGS brightness, as measured by the flux incident on the AO wavefront sensor, 
was equivalent to a 1/ ~ 9.2 — 10.3 mag star. Tip-tilt aberrations and quasi-static changes in the 
image of the LGS as seen by the wavefront sensor were measured contemporaneously with a second, 
lower-bandwidth w avefront sensor monitoring the R = 16.2 mag field star USNO-Bl.O 0601-0344964 
(JMonet et ahlboosl l. located 31" away from 2MASS J1534-2952AB. 



At each epoch, 2MASS J1534— 2952AB was imaged in filters covering the sta ndard 2.2 um at- 
mospheric window from the the Mauna Kea Observatories (MKO) filter consortium (ISimons &: Tokunaga 
2002; iTokunaga et al.l 12002) . Our initial observations in April 2005 were carried out with the K' 
(2.12 //m) filter to minimize the thermal background from the AO system, which is kept at ambient 
temperature. Subsequent runs employed the K (2.20 //m) or Ks (2.15 /xm) filters. Hereinafter, for 
brevity we refer to all these data simply as i^-band observations. 

On each observing run, we typically obtained a series of dithered i('-band images, offsetting 
the telescope by a few arcseconds between each 1-2 images. The sodium laser beam was pointed at 
the center of the NIRC2 field-of-view for all observations. In April 2005, we also obtained images 
with the MKO J (1.25 /xm) and H (1.64 /im) filters. In April 2008, we also obtained images with 
the CH4S filter, which has a central wavelength of 1.592 fim and a width of 0.126 /xm; this filter is 
positioned around the iJ-band fiux peak in the spectra of mid/late-T dwarfs. 



The images were reduced in a standard fashion. We constructed flat fields from the differences 
of images of the telescope dome interior with and without continuum lamp illumination. Then 
we created a master sky frame from the median average of the bias-subtracted, flat-fielded images 
and subtracted it from the individual images. Images were registered and stacked to form a final 
mosaic, though all the results described here were based on analysis of the individual images. 
Outlier images with much poorer FWHM and/or Strehl ratios were excluded from the analysis. 
Instrumental optical distortion was corrected based on analysis by B. Cameron (priv. comm.) of 
images of a precisely machined pinhole grid located at the first focal plane of NIRC2. The la 
residuals of the pinhole images after applying this distortion correction are at the 0.6 mas level over 
the detector field of view. Since the binary separation and the imaging dither steps are small, the 
effect of the distortion correction is minor, smaller than our final measurement errors. 

Table [T] compiles the details of our observations, and Figure [T] presents our Keck LGS data. 
Full widths at half maxima (FWHM) and Strehl ratios were determined from two field stars located 
f«5-6" from 2MASS J1534-2952AB. The tabulated errors on the FWHM and Strehl ratios are the 
standard deviation of these quantities as measured from the individual images. 

To measure the flux ratios and relative positions of 2MASS J1534— 2952AB's two compo- 
nents, we mostly used the two aforementioned nearby field stars, observed simultaneously with 
2MASS J1534— 2952AB on NIRC2. These stars provided an excellent measurement of the instan- 



taneo us PSF. We empirically modeled the PSF using the Starfinder software package (JDiolaiti et al. 



2000), which is designed for analysis of blended AO images. For the Jan 2008 data, we employed 
a different procedure, fitting analytic PSFs comprising multiple elliptical gaussians to model the 
binary images. These data were taken at much higher airmass than all the other data. Be- 
cause of the larger atmospheric dispersion and the different colors of the field stars relative to 
2MASS J1534— 2952AB, PSF fitting produced less accurate results than the analytic approach, as 
determined by the artificial binary tests described below. For every image, we fitted for the fiuxes 
and positions of the two components and then computed the flux ratio, separation, and position 
angle (PA) of the binary. The averages of the results were adopted as the flnal measurements. 
Overall, our PSF fitting produced very high quality relative measurements, with errors of order 1% 
for the flux ratios, 0.05 pixels for the binary separation, and 0.2° for the PA. Note that these latter 
two values account only for the internal instrumental measurements and do not include the errors 
on the astrometric calibration of NIRC2, which we include below. 

In order to gauge the accuracy of our measurements, we created myriad artificial binary stars 
from images of the two PSF stars. One PSF star was used to create artificial binaries, and the 
other was used as the single PSF for fitting the components. For data at each epoch, Starfinder was 
applied to the artificial binaries with similar separations and fiux ratios as 2MASS J1534— 2952AB. 
These simulations showed that any systematic offsets in our fitting code are very small, well below 
the random errors, and that the random errors are accurate. In cases where the RMS measure- 
ment errors from the artificial binaries were larger than those from the 2MASS J1534— 2952AB 



measurements, we conservatively adopted the larger errorso 



To convert the instrumental measurements of the binary separation and PA into celestial 
units, we used a weighted average of the calibration from lPravdo et alj (J2006l ). with a pixel scale of 
9.963 lb 0.011 mas/pixel and an orientation for the detector's +y axis of —0.13 it 0.07° east of north. 
These values agree well with Keck Observatory's notional calibration of 9.942 ±0. 05 mas/pixel and 



0.0 ±0 .5°, as well as the 9.961 ±0.007 mas/pixel and -0.015 ±0.134° reported bv lKonopackv et al 



(|2007l ). Also, comparison of NIRC2 images of M92 to astrometrically calibrated HST/ ACS Wide- 
Field Camera images gives a pixel scale for NIRC2 that agrees to better than 1 part in 10"'^ with 
our values (J. Anderson, priv. comm.). 

Finally, we must consider the effect of atmospheric refraction. Because of the southern dec- 
lination of 2MASS J1534— 2952AB, all of our Keck observations were necessarily undertaken at 
significant airmass (>1.55). Because the two components of the binary do not have exactly the 
same spectral types (§ 13. ip . the observed positions on the sky are subject to slightly different 
amounts of differential chromatic refraction (DCR). We computed t he exp ected shift in the relative 



astrometry at each epoch using the prescriptions of iMonet et al.l (|l992l ) for the DCR offset and 



Stond (| 19841 ) for the refractive index of dry air. We assui ned a fiducial teniperat ure of 275 K and 
pressure of 608 millibars for conditions on Mauna Kea (ICohen &: Cromerill988l ). We computed 
the effective wavelengths for spectral types of T5.0 and T5.5 for the two components using us- 
ing all available spec t ra of these subc lasses contained in the SpeX Prism Spectral Library (from 
Burgasser et al.ll2004l : IChiu et al.ll2006l . and lLooper et al.ll2007l ) and the appropriate filter response 
curvecl Note that because of the unusual spectra of T dwarfs, the effective wavelengths of the 
secondary component is actually bluer in H- and i^-bands and redder only in J-band compared 
to the primary component. Given the fact that the secondary is mostly north of the primary at 
all of our Keck epochs, DCR causes the separation of the binary to appear slightly smaller at H 
and i^-bands and slightly larger at J-band compared to the true position as would be observed at 
zenith. The amplitude of the DCR effect is about 0.3 mas, much smaller than the measurement 
errors at most (but not all) of the Keck epochs. However, the effect is a systematic one so we 
correct the relative astrometry of the two components based on our calculations. 



■^We also experimented with directly fittin g the binary data by itself, without us ing other stars as the PSF. This 
approach was similar to our previous analyses ( Liu fc Leggettll2005l : iLiu et al.ll2006a ). namely we fit the images of the 
binary with either (1) an analytic model of the point spread function (PSF) as the sum of elliptical Gaussians, or (2) 
an empirical model derived iteratively using Starfinder-based code. These were also tested against images of simulated 
binaries. As before, we found that the Starfinder measurements were somewhat better compared to the analytical 
fits for the iiT-band data, and the multi- Gaussian fitting was better for the J and H-band data. Overall, fitting the 
binary images by themselves produced very good astrometry — with internal errors (i.e., without the uncertainty in 
the absolute NIRC2 astrometric calibration) of about 3% for the flux ratio, 0.1 pixels for the separation, and 0.3° 
for the PA — but slightly worse than our fits using field stars as PSFs. The exception was the Jan 2008 data, as 
described in the text. 



MKO filter curves available at ftp://ftp.jach. hawaiLedu/pub/ukirt/'skl/filters 



Table 1 presents the final resulting measurements from our Keck LGS data. For the April 2005 
dataset, all three filters give astrometry consistent within the measurement errors; we use only the 
iif-band results in the orbit fitting discussed below. In the Table and in our orbit fitting (§ 13. 3p . 
we take care to discriminate between the instrumental errors (namely those that arise solely from 
fitting the binary images) and the overall astrometric calibration of NIRC2, and thus any future 
refinements in the latter can be readily applied to our measurements. 



2.2. HST 

2.2.1. WFPC2 Planetary Camera 

The two components of 2MASS J1534-2952AB are only barely resolved in the ifST/WFPC2 
F^lAW discovery images from August 2000. Therefore, to determine their relative positions and 
fluxes, we must model the images using the sum of two blended PSFs. The PSF of WFPC2's 
Planetary Camera (PC) is undersampled (FWHM = 1.7 pix for FSl^W); this makes any empirical 
determination of the PSF difficult without PSFs sampled at many subpixel locations. Moreover, 
Anderson & King (2003) found that the WFPC2 PSF varies significantly over the detector due 
to geometric distortion, making it impossible to construct a reliable empirical PSF from other 
stars in th e same image, ev e n if the re are enough to sample many subpixel locations. The original 



analysis by lBurgasser et al.l ()2003bl ) employed a hybrid gaussian/empirical PSF to fit for the binary 
parameters with resulting uncertainties of ±7 mas in separation and ±9° in PA. The astrometry 
from the WFPC2 discovery epoch is obviously very important to the orbit determination, so we 
undertook our own analysis with a more accurate PSF model to improve the precision of the binary 
parameters. 

We used the TinyTim software package (Krist 1995) to create model PSFs for the WFPC2 
images. We generated 5x super-sampled PSFs that included the effects of (1) variation with 
position on the detector; (2) broad-band wavelength dependence, by taking into account the filter 
response curve and the spectrum of the source (using the Keck/LRIS optical spectrum of the T4.5 
dwarf 2MASS J05591914- 1404488 from Burgasser et al. (2003) as the template for the individual 
components of 2MASS J1534— 2952AB), (3) telescope jitter (0 to 20 mas of gaussian jitter), and 
(4) telescope defocus (±10 ^m) to account for HST breathing effects. Because the geometrical 
distortion is location-dependent, we used TinyTim model PSFs generated for the nearest integer 
pixel location to the centroid of the binary or single T dwarf. Also, we used the template spectrum 
closest to the spectral type of the T dwarf with sufficient wavelength coverage (0.70-0.96 jjlth) from 
S. K. Leggett's spectral library^ 



^The optical spectra we used were: 2MASS J05591914- 1404488 fo r T4.5 and T5.5 objects (JBurgasser et al.ll2003ah : 



SDSSp J162414.37+00 2915.6 for T6 objects (JBurgasser et al.ll200ol): SDSSp J134646.45-003150.4 for T6.5 objects 



( Burgasser et al.ll200ol ): and Gl 570D for T7.5 objects (JBurgasser et al.ll2003al l 



These TinyTim model PSFs were used to fit simultaneously for (1) the location of the primary, 
(2) the location of the secondary, (3) the normalization of the model PSF to the primary, and 
(4) the flux ratio of the two components. When fitting positions, the super-sampled TinyTim 
PSF was interpolated using cubic convolution to the appropriate subp ixel location. The best fit 



values were found using the amoeba alogorithm (e.g., iPress et al.l 11992 ) to find the minimum x^ 
value of a 1.1" x 1.1" subimage centered on the binary. The image was cleaned using the IDL 
routine CR_REJECT in the Goddard IDL library to identify and mask the numerous cosmic rays in 
the undithered WFPC2 image pair. Masked pixels were excluded from the computation of the 
X^ value. The noise in each pixel was determined from the bias-subtracted raw WFPC2 image, 
assuming a read noise of 5.3 e~/pix and a gain given by the header keyword ATODGAIN. A grid of 
PSFs in telescope jitter and defocus were tried, and the fit corresponding to the jitter and defocus 
combination yielding the lowest x^ was chosen. For images of 2MASS J1534— 2952AB we found 
that our PSF-fitting routine yielded residuals of <2% of the peak value in 90% of pixels with 
S/N > 3. 

Due to the optical distortion and "34-th row" defect present in the WFPC2 (Anderson &; King 
1999, 2003), the best-fit pixel locations of each binary component do not exactly correspond to 
their locations in an undistorted celestial reference frame. To remove these effects, we applied the 
corrections of Anderson &; King (2003), using a pixel scale of 45.54 it 0.01 mas/pix|£| The 34-th row 
effect could change, for example, the binary separation by as much as 0.7 mas (a systematic shift 
of — 0.6(t), but because the binary components do not straddle a defective row and are separated 
by a mere 1.4 pixels, the application of the 34-th row and distortion corrections have a negligible 
effect on the astrometry. 

With only two undithered WFPC2 images of 2MASS J1534-2952AB, it is challenging to 
quantify the measurement uncertainties, and it is impossible to completely quantify the systematic 
errors, which arise from an imperfect PSF model and also probably depend on the subpixel positions 
of the two components given the undersampled nature of the data. Using only the RMS scatter of 
the two measurements, the inferred random errors in separation, PA, and flux ratio are 0.9 mas, 
0.07°, and 0.03 mag, respectively. To derive more robust random errors and to investigate the 
systematic errors, we conducted an extensive Monte Carlo simulation of our fitting routine. 

We used WFPC2 FSIAW images of seven other T dwarfs from the same HST program, 
all apparently single, to create artificial binaries which we then modeled using our PSF-fitting 



^The value of the pixel scale and its uncertainty come from the WFPC2 Instrument Handbook for Cycle 13. This 
number is consistent with other measurements available in the literature: (1) Holtzman et al. (1995) derived a pixel 
scale of 45.54 mas/pix by comparing commanded telescope offsets in arcsec to the resulting pixel offsets; (2) Holtzman 
et al. also derived a pixel scale of 45.55 mas/pix by comparison to an astrometric standard field in M67; (3) Pascu et 
al. (1998) used the JPL ephemeris of the satellites of Uranus to derive a pixel scale of 45.57 mas/pix; and (4) Pascu 
et al. (2004) used the JPL ephemeris of the satellites of Neptune to derive a pixel scale of 45.55 mas/pix. The scatter 
in these pixel scales is consistent with our quoted uncertainty. 



-9- 

routinejil Because the WFPC2 PSF is severely undersampled, we only created artificial binaries 
with integer-shifted positions. It turns out that the WFPC2 locations of the two components of 
2MASS J1534— 2952AB are at a very nearly integer-shifted separation of 1.4 pixels (Ax ~ 1 pix, 
Ay ~ — 1 pix). Therefore, in determining the uncertainties and systematic offsets we used only the 
configuration most nearly matching that of 2MASS J1534— 2952AB with a separation of \/2 pix 
and instrumental PA of 225°. We found that using any or all of the other three v2 configurations 
gave consistent uncertainties. After subtracting the best-fit model, artificial binary images yielded 
residual images in which 90% of pixels with S/N > 3 were below 1-3% of the peak flux, comparable 
to the residual images of 2MASS J1534-2952AB. 

The images of the single T dwarfs span a range in S/N from about 1.5 mag brighter to 
1.3 mag fainter than the primary component of 2MASS J1534— 2952AB. We used these images at 
their native S/N when simulating the primary component. To simulate the secondary component, 
we degraded the S/N of the images assuming a flux ratio of 0.30 mag. We also tried flux ratios 
of 0.25 mag and 0.35 mag to explore the possibility that the uncertainties depend on the assumed 
flux ratio, but we found that this had an insignificant effect on our predicted uncertainties (< la^, 
where a^ = a/y/2Nsim and Ngim is the number of simulations). Signal-to- noise degradation of an 
image was done by a multiplicative scaling followed by the addition of normally distributed random 
noise to each pixel, according to the same WFPC2 noise model we used to determine x^ ^^ the PSF- 
fitting procedure. In fact, by running simulations where the primary images were degraded to much 
lower S/N, we found that all of the single T dwarfs are in a high S/N regime in which systematic 
errors (PSF model imperfections) dominate, not random errors (photon noise) — our simulations 
showed no dependence between the S/N of the T dwarf used to construct artificial binaries and the 
resulting astrometric uncertainties. Therefore, we used the RMS of the results from all simulated 
binaries in order to determine the final uncertainties for 2MASS J1534— 2952AB. As expected, the 
uncertainties in separation, PA, and flux ratio from our simulations were somewhat larger than those 
derived from the standard deviation of the measurements from the two 2MASS J1534— 2952AB 
images, since both random and systematic errors have been evaluated in the simulations. In fact, 
the larger uncertainties are not simply due to averaging over, for example, the many subpixel 
locations of the single T dwarfs used in the Monte Carlo, because the simulated measurements for 
each single T dwarf show scatter consistent with the final derived uncertainties. 

Table [2] presents our final results for the WFPC2 images, with systematic offsets from the 
Monte C arlo simulations a pplied . Our astrometry agrees well (better than lo") with the original 



results of lBurgasser et al.l (|2003bl ). though our measurement errors are a factor of 8 smaller. Part 
of this improvement comes from our use of TinyTim-computed PSFs, as opposed to the simpler 
Burgasser et al. PSF model of a gaussian plus empirical residuals. We also used all possible 



■^These were 2MASS J05591914-1404488 (T4.5), 2MASSI J0937347+293142 (T6.0), 2MASSI J1217110-031113 
(T7.5), 2MASS J12373919+6526148 (T6.5), Gl 570D (T7.5), 2MASSI J15 46291-332511 (T5.5), and 
2MASSI J2356547— 155310 (T5.5), using spectral types from iBurgasser et al.l ( 2006a ). The remaining object from 



tliis program, 2MASSI J1047538+212423 (T6.5), was unusable for our purposes because it landed on a bad column. 
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single PSFs in our artificial-binary simulations, whereas Burgasser et al. used only WFPC2 images 
of 2MASS J0559— 1404, a source that is suspected to be an unresolved binary|f| And part of 
the improved uncertainties is somewhat illusional, as it arises from the different parameter space 
explored in Monte Carlo simulations by us and Burgaser et al. As a check, we ran a suite of 
simulations more comparable to that of Burgasser et al., in which the ranges of artificial binary 
parameters were 1.0-3.0 pixels in separation, the full range of PAs, and 0.0-1.0 mag in flux ratio. 
These yielded similar astrometric uncertainties to the published errors, suggesting that the smaller 
uncertainties we derive are due to our more restricted choice of artificial binary configurations 
(Ax = 1 pix. Ay = — 1 pix) and/or averaging over many single PSFs to reduce systematics 
associated with any one specific object. 

Our improvement to the WFPC2 astrometry was essential in our early attempts to fit the orbit 
based on Keck data obtained in 2005-2007. However, with the addition of data in 2008, the final 
ifST+Keck dataset has sufficent time baseline and astrometric quality that the choice of WFPC2 
astrometry does not highly impact the final orbit fitting results (§ l3.3.3p . 



2.2.2. ACS High Resolution Camera 

The 2MASS J1534-2952AB system was observed on 2006 January 19 and 2006 April 11 (UT) 
with the High Resolution Camera (HRC) of HSTs Advanced Camera for Surveys (ACS) by program 
GO-10559 (PI H. Bouy). The binary is much more widely separated at these epochs than in the 
WFPC2 observations, but the PSFs of the two components are still blended. We have therefore 
applied the same TinyTim PSF-fitting technique described in the previous section to derive the 
relative astrometry from the ACS images. The primary differences between the WFPC2 and ACS 
datasets are: (1) ACS has much more severe geometric distortion than WFPC2, which changes 
the shape of the PSF and complicates astrometry because the pixels projected on the sky are not 
square, and (2) the ACS data are of much lower S/N ., with a total exposure time of only 50 sec 
for each cosmic-ray rejected, combined dithered image (cf., 1300 sec for a single WFPC2 image). 
Because of the lower S/N , we found it unwarranted to fit the ACS images for telescope jitter and 
defocus as adding these free parameters did not improve the quality of the fits (as verified in the 
Monte Carlo simulations discussed below). Also, we found that the ~25x lower S/N of these data 



^2MASS J0559— 1404 is rquRhly twice as luminous as objects of similar spectral type I Dahn et al.1 12002 : 



Tinne^_et_^ 2003 : Vrba et al.M2004l l. High angular resolution observations have not detected any multiplicity 
IjBurgasser et al.ll2003bl : Liu et al., in prep.), but the source could be a very tight system, e.g., with a 0.5 pixel 



separation. If so, its multiplicity could confuse any attempts to fit only two single PSFs to the tightest artificial 
binaries constructed from its image. Interestingly, we found that at the smallest (1.0 pixel) separations, artificial 
binaries constructed from images of 2MASS J0559— 1404 yielded extremely large uncertainties (5 mas and 15°). We 
did not observe such behavior for any of the other six apparently single T dwarfs, nor were the uncertainties for 
larger separation binaries made from 2MASS J0559— 1404 abnormally large. One natural explanation would be that 
the source is just marginally resolved in the WFPC2 imaging. 
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almost exactly negates any improvement to the astrometry that might be expected given the larger 
binary separation at these epochs. 

We used distorted model PSFs generated by TinyTim to fit for the position and flux of each 
binary component in images that had been cosmic-ray cleaned (CRSPLIT=4) by the latest HST 
pipeline. Best-fit pixel locations were corrected for geometric distortion using the solution of An- 
derson &; King (2004; Instrument Science Report 04-15), and we used their measured ACS pixel 
scale, which was derived by comparing commanded (POSTARG) offsets of HST in arcseconds to the 
resulting pixel offsets. They derived two such pixel scales for two epochs of observations of 47 Tuc, 
and we adopt the mean and standard deviation of these two: 28.273 it 0.006 mas/pix|£| For each 
epoch, we adopted the mean of the measurements from all four dithered images for the binary 
parameters of 2MASS J1534-2952AB (Table E]). 

Again, to investigate the measurement errors thoroughly, we performed Monte Carlo simula- 
tions of our fitting routine. We used images of single brown dwarfs to construct artificial binaries in 
configurations resembling 2MASS J1534— 2952AB. At both ACS epochs, the binary is well repre- 
sented by integer-pixel shifts on a grid where Ax = (6, 7) and Ay = (—3, —2, —1,0, 1). There have 
been no HST/ ACS science programs dedicated to studying single b rown dwarfs; however, ACS 



images of the single brown dwarfs 2 MASS J00361617+ 182 1 104 (L3.5; iKirkpatrick et al.ll2000l ) and 



2MASS J05591914- 1404488 (T4.5: iBurgasser et al.ll2006bl ) were obtained for calibration purposes 



and are available in the IfST Archive (CAL/ACS-10374, PI Giavalisco). We found that despite the 
large difference in spectral types, any corresponding difference in the PSF does not alter the results 
of the Monte Carlo simulations. The S/N of each of these single objects is much higher than that 
of 2MASS J1534— 2952AB, so we degraded the S/N of the single brown dwarfs for the artificial 
binary simulations. In fact, by varying the S/N of the simulations, we found that the ACS data for 
2MASS J1534— 2952AB are well in the S/N regime dominated by random photon noise, while the 
images of the single objects are in a high S/N regime dominated by systematic errors (akin to the 
WFPC2 images of 2MASS J1534— 2952AB). Therefore, given that we have four images, we divide 
the RMS of the Monte Carlo results by v4 to represent the final uncertainties. 

Table [2] contains our final ACS results!^ Note that the ACS data are contemporaneous with 
our Keck LGS data, and observations by the two telescopes in 2006 separated by less than one 
month show excellent agreement. However, the relatively low S/N of the ACS data means it has 
the larger astrometric errors. Our simulations confirm that if the ACS data were of higher S/N 
(i.e., longer exposure times than 50 s), the resulting astrometric precision would be much better 



^Comparison of WFPC2-PC images of 47 Tuc (GO-8267, PI Gilliland) with ACS-WFC images (GO-10775, PI 
A. Sarajedini) shows excellent agreement between our adopted pixel scales for the two instruments, at the level of 
2 X 10^'* which is well below the other errors in the measurements (J. Anderson, priv. comm.). 

^ After our paper was submitted, IBouv et al.l ( 2008 ) reported an analysis of the same ACS images. Their results 



10, 



agree with ours to within the stated uncertainties. Their errors are slightly smaller than ours in separation (0.9 mas 
compared to 1.1 mas) and much smaller in PA (0.1° compared to 0.5°). 
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than, instead of comparable to, the WFPC2 astrometric uncertainties. 



2.3. IRTF/SpeX Spectroscopy 



We obtained low-resolution {R ?»150) integrated-hght spectra of 2MASS J1534-2952AB on 
2008 May 16 UT from NASA's Infrared Telescope Facility (IRTF) located on Mauna Kea, Hawaii. 
Conditions were photometri c with seeing of abo ut 0.7" FWHM near the target. We used the facility 
near-IR spectrograph Spex (JRavner et al.lll998l ) in prism mode, obtaining 0.8-2.5 ^m spectra in a 
single order. We used the 0.5" wide slit, oriented at the parallactic angle to minimize the effect of 
atmospheric dispersion. 2MASS J1534— 2952AB was nodded along the slit in an ABBA pattern, 
with individual exposure times of 180 sec, and observed over an airmass range of 1.64-1.60 as 
it rose. The telescope was guided during the exposures using images obtained with the near-IR 
slit-viewing camera. The total on-source exposure time was 720 sec. We observed the AO V star 
HD 142851 contemporaneously for flux and telluric calibration. All spectra were reduced using 
version 3.4 of the SpeXtool software package (IVacca et al.ll2003l : ICushing et al.ll2004l ). The reduced 
IRTF /Spex spe c trum i s plotted in Figure [2] and compared to T dwarf spectral standards from 
Burgasser et al.l (l2006bl ). Visual examination shows an excellent match to the T5 spectral standard 
2MASS J1503-I-2525, as does measurement of the iBurgasser et al.l (|2006bl ) spectral indicies for 
2MASS J1534-2952AB: H2O-J = 0.271 (T4.8), CH4-J = 0.420 (T4.8), R2O-H = 0.345 (T5.0), 
CRi-H = 0.430 (T5.0), and CR^ -K = 0.22 4 (T5.1 '), with spectral type estimates based on the 



polynomial fits to the indices from [BurgasserJ (|2007al ). 



3. Results 



3.1. Resolved Photometry and Spectral Types 



We use our measured flux ratios and the published JHK photometry from lKnapp et al.l (120041 ) 
to derive resolved IR colors ar id magnitudes for 2MAS S J1534— 2952AB on the MKO system. We 
use the HST photometry from IBurgasser et al.l (|2003bl ) in determining the resolved FSIAW magni- 
tudes. Then to infer spectral types for the individual c omp onents, we compar e these to magnitudes 



and colors of ultracool dwarfs from lKnapp et al.l (|2004l) and 



Chiu et al 



20061) . excluding known bi- 



naries. We use near-IR spectral classfications from the IBurgasser et al.l (J2006bl ) scheme. We assume 
that the components of 2MASS J1534— 2952AB are themselves single, not unresolved binaries. 

Figure [3] shows that component A has IR colors most typical of T4.5-T5 dwarfs, and compo- 
nent B is most similar to T5-T6 dwarfs. The individual absolute IR magnitudes (given in Table [3]) 
give similar res u lts. Th e "faint" polynomial fits for absolute magnitude as a function of spectral type 
from lLiu et Zl (|2006bl ) give M{J) = {14.5, 14.6, 14.7, 14.9, 15.2} mag, M{H) = {14.6, 14.8, 15.0, 
15.2, 15.5} mag, and M{K) = {14.7, 14.9, 15.1, 15.4, 15.7} mag for near-IR spectral types of T4.5, 
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T5, T5.5, T6, and T6.5, respectively. Averaging the same data for each individual subclass gives 
M{J) = {13.9±0.6, 14.1, 14.4±0.4, 15.0±0.5, 15.1±0.5} mag, M{H) = {14.0±0.6, 14.2, 14.6±0.4, 
15.3±0.4, 15.4±0.5} mag, and M{K) = {14.0±0.5, 14.3, 14.6±0.4, 15.5±0.7, 15.6±0.9} mag, where 
the uncertainties are the RMS of the photometry for each subclass (and no listed uncertainties for 
subclasses with only one object). Altogether, the absolute magnitudes suggest types of T5-T5.5 
for component A and T5.5-T6 for component B. 

The resolved {CH/j^s — H) colors provides a third means to estimate the spectral types, as 
these track the H-hand me thane absorpt i on, w hich correlates well with overall near-IR spec- 
tral type (e.g.. Figure 2 of iTinnev et al.l |2005| ). First, we compute individual CH^s magni- 
tudes for 2MASS J1534— 2952AB using the CH^s fl ux ratio from ourLGS images, the integrated- 
light photometry oi H = 14.74 it 0.03 mag from iKnapp et al.l ()2004l ). and an integrated- light 
color o f {CH/is — H) = —0.33 mag synthesized from the near-IR spectrum of iBurgasser et al. 



(J2006bl ). Including the measurement errors in the flux ratios and H-hand photometry, we find 
CH^s = 15.06 ± 0.04 and 15.27 ± 0.04 mag and [CHiS - H) = -0.30 ± 0.05 and -0.37 ± 0.05 mag 
for components A and B, respectively. Note that the relative {CH^s — H) color of the two compo- 
nents is known to higher precision, since the above computed colors for the two components contain 
the same 0.03 mag error that originates from the integrated-light H-hand photometry. (In other 
words, the 0.05 mag uncertainties in the colors of the two components are not independent, but 
correlated.) Removing this effect gives a relative color of A.[CHiS — H) = 0.07 it 0.02 mag between 
A and B, i.e., greater methane absorption in component B is detected. 

To determine the behavior of [CHiS—H) with near-IR spectral type, we synthesized colors from 
the Spex Prism Spectral Library collection, which contains low-resolution spectra of 68 T dwarfs 
after removing spectrally peculiar objects and known binaries^ Figure S] plots the results. We fit 
a 2nd-order polynomial for the dependence of color on near-IR spectral type and vice-versa: 



{CHiS -H) = -0.400 + 0.07718 x SpT - 0.0029736 x SpT'^ 

SpT = 19.40 - 19.698 x {CH^s - H) - 8.3600 x {CHiS - Hf 



(1) 
(2) 



where SpT = 20 for TO, = 21 for Tl, etc. The RMS scatter about the fits are 0.02 mag and 
0.3 subclasses, respectively. Using these polynomial relations, the observed (CH^s — H) colors give 
spectral types of T4.5 it 0.7 and T5.6 it 0.6 for components A and B, respectively, where the spectral 
type uncertainties come from formal propagation of the uncertainty in the colors. In addition, just 
as the relative {CH4S — H) color of the two components are known more accurately than the 
absolute colors, we compute a relative spectral type of 1.1 it 0.3 subclasses between components A 
and B. 

Combining the inferences from the JHK colors, the absolute magnitudes, and the {CH4S — H) 
colors, we adopt spectral type estimates of T5 it 0.5 and T5.5 it 0.5 for the two components. The 



jhttp : //www . br owndwarf s ■ org/spexprism 



20061 



nttp: //www. browndwaris.org/spexprism . 
bllali: lBurgasseiJ l|2007bl ): IChiu et all ( |2006l ) 



The T dwa r f spe ct ra are compiled from Burgasser et al.1 l|2004 
: iLooper et all (|2007l ') : iLiebert fc Bmgassej (|2007l '). 
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relative {CHAg — H) color favors a slightly larger spectral type difference than the absolute magni- 
tudes but consistent with the adopted uncertainty. (Also the T^n difference of the two components 
computed in ^ 14.31 favors a 0.5 subclass difference.) 



Higher-orde r multiple systems are very rare among ultracool binaries, with an estimated fre- 
quency of 3l^% (JBurgasser et al.ll2006dl ). and thus a priori we do not expect 2MASS J1534— 2952AB 
to fall into this category. The colors and magnitudes are consistent with the system being com- 
posed of only two components, and not being a partially resolved higher order multiple system. If 
component B was actually an equal-mass binary, the absolute magnitudes of its components would 
be 0.75 mag fainter than the integrated-light of B, meaning Ad{J, H, K) = {15.5, 15.9, 16 .0) mags 



This w ould suggest a spectral type around T7, based on the polynomial relations in iLiu et al 



(|2006bl ). which is clearly too late- type compared to the integrated-light spectrum and the observed 
near-IR colors of B. 



3.2. Bolometric Luminosities 



To measure the Lboi for the system, we combine our SpeX 0.9-2.4 /im spectrum with the 
publishe d integrated-light F814VF, L^-band, and Spitzer/IRAC thermal-IR photometry and unc er- 



tainties (JBurgasser et al.ll2003bl : lKnapp et al.ll2004l : iGolimowski et al. 



2004a 



Patten et al. 



200' 



We flux-calibrated the SpeX data using the published H-haxid MKO photometry from Knapp et 
al. 2004. For the Spiteer/IRAC data, we adopted the photometric calibration and an overall 2% 
absolute uncertainty based on Reach et al. (2005). We extended the binary's spectral energy dis- 
tribution (SED) to short wavelengths by linearly interpolating from the FSIAW datum to zero flux 
at zero wavelength and to long wavelengths by assuming a Rayleigh-Jeans spectrum beyond the 
reddest Spitzer bandpass (7.87 /xm); this extrapolation increases the total flux by 2%. 

We then integrated the SED, using a Monte Carlo approach to account for all the measurement 
errors. We find log(Lboi/-^0) = — 4.751ib0.011 dex for the system, with the uncertainty increasing to 
0.018 dex after including the uncertainty in the distance. (As discussed in § 4, we keep track of these 
independent uncertainties in our calculations.) The largest uncertainty in the integration arises from 
the 0.03 mag uncertainty in the integrated-light photometry used to normalize the SpeX spectrum. 
We cross-checked our method using the same data f or the T4.5 dwarf 2M ASS J0559— 1404 and 
found excellent agreement with the Lboi measured by ICushing et al.l (120061 ) using absolutely flux- 
calibrated spectra from 0.6-15 /im. 



^■^The L'-band photometry for this system is anomalous, as iLeggett et al.1 ( 2007a ) show the ([3.5] — L') color is 
~0.2 mag redder than any other T dwarf and ~0.3 mag redder than objects of similar spectral type. Both of these 
bandpasses are affected by the 3.3 /im fundamental band of CH4, with the Spitzer [3.5] data being subjected to greater 
absorption, and thus the redder color might point to anomolously strong CH4. However, the near-IR spectrum and 
the JHK and Spitzer 3.5-7.9 /xm colors are consistent with other mid-T dwarfs, and thus perhaps the L'-band 
photometry is incorrect. Either way, this has a negligible effect on the computed Lboi- 
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The computed total Lboi agrees well with that inferred from using the i^-band bolometric 
corrections [BCk) from Golimowski et al. (2004), namely using the resolved iiT-band absolute 
magnitudes and the estimated spectral types, which would give log(Lboi/-^0) = —4.97 it 0.06 and 
—5.06 lb 0.06 dex for the individual components and thus log(I/boi/-^0) = —4.71 ± 0.08 dex for 
the total system. However, the uncertainties are larger when using BCk to derive -Lbob since this 
incorporates the uncertainties arising from the 0.5 subclass uncertainty (0.06 mag in bolometric 
magnitude) and the intrinsic scatter in the Golimowski et al. BCk relation (0.13 mag). In short, 
direct integration of the observed SED is more accurate. 

To apportion the observed total Lboi into the individual components, we assume the observed 
iiT-band flux ratio of the system represents the luminosity ratio. This would be exactly correct 
if the two components had identical spectral types (and neglecting photometric variability). To 
account for the difference in spectral types, we generate a Monte Carlo distribution of BCk values 
for each component using the Golimowski et al. polynomial fit as a function of spectral type subject 
to the following rules: the spectral type of component A is uniformly distributed from T4.5-T5.5; 
the spectral type of component B is no later than T6; and the difference in their spectral types is 
at least 0.5 subclasses. This produces an average difference in BCk between the two components 
of 0.09 mags and an RMS of 0.03 mags. Thus, we find log(Lboi/^0) = -5.015 ± 0.019 and 
—5.093 lb 0.019 dex for the two components, including the uncertainty in the distance. 



3.3. Dynamical Mass Determination 

3.3.1. Orbit Fitting using Markov Chain Monte Carlo 

We have data at 9 independent epochs, which is formally sufficient to determine the 7 pa- 
rameters of a visual binary orbit given our measurements (9 positions + 9 times). However, two 
pairs of measurements are separated by only one month (April/May 2006 and March/April 2007) , 
and the cadence of the orbital phase covered is limited, with the iJST/WFPC2 datum being taken 
almost 5 years before the next epoch. While standard gradient-descent (Levenberg-Marquardt) 
least-squares techniques would be sufficient to derive an orbit (^ I3.3.2p . we also would like to accu- 
rately determine the probability distribution of the orbit parameters (which may not be normally 
distributed) and the associated degeneracies. For epochs with Keck data taken in multiple filters, 
we choose the filter with the smallest astrometry errors. 

Thus, we first used a combination of gradient-descent techniques from random starting points 
and simulated-annealing techniques to isolate the class of potential orbital solutions near a reduced 
chi-squared (x^) of 1. Then, to fully explore thi s class of solutions, we used a Markov Chain Monte 



Carlo (MCMC) approach (e.g. iBremaudI Il999l ). MCMC provides a means to explore the multi- 
dimensional parameter space inherent in fitting visual orbits that is computationally efficient, able 
to discern the degeneracies and non-gaussian uncertainites in the fit, and allows for incorporation of 
a priori knowledge. In short, the MCMC approach is distinct from ordinary Monte Carlo methods 
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in that instead of a completely random steps through the model parameter space, the steps are 
chosen such that the resulting number of samples (the "chain" ) is asym ptotically equivalent t o the 



poste rior probability dis tribution of the parameters being sought. (See lTegmark et al.ll2004l . iFord 



20051 . and lGregorvll2005l for explications of applying MCMC to astronomical data. 



We parameterized the binary's orbit using the standard 7 parameters: period (P), semi-major 
axis (a), inclination (i), epoch of periastron (Tq), PA of the ascending node (Ori. eccentricity 
(e), and the argument of periastron (lo). We used an MCMC chain length of 2 x 10^, with the 
parameters stored every hundredth iteration. We started the chain at the global minimum found 
by the simulated-annealing and gradient-descent algorithms. We used the Metropolis-Hastings 
algorithm to sample the joint probability distribution with a variant of the usual Gibbs sampler. 
Instead of chosing one of the 7 parameters to increment or decrement, we chose instead to move 
randomly either forwards or backwards along one of 7 orthogonal directions in the parameter 
space. These directions were initially chosen to be along the individual coordinate axes of our 
7 parameters. Then, every 5 x 10^ iterations, the covariance matrix of the 7 parameters was 
calculated and new directions chosen along the directions corresponding to the eigenvectors of the 
covariance matrix. After covariance matrix calculation, the trial step size was set to be the square 
root of the covariance matrix eigenvalues. This enabled the long, thin, curved minima in our 
parameter space to be sampled much more efficiently. Each set of 5 x 10^ iterations can thus be 
thought of as its own sub-chain with fixed directions. As the first 5 x 10^ iterations uses a less 
efficient sampling (sampling directions aligned with the 7 parameters themselves), we treated this 
part of the chain as the "burn-in" time, neglecting it in the final analysis. As the chain was running, 
the trial step size in each direction was continuously scaled over a timescale of 700 steps so that 
the success r ate of jumps ayeraged to 0.25. The correlation length of our most correlated chain. 



as defined in iTegmark et al.l l2004l . was 2 x 10^ for the orbital period (P), with equal or smaller 
correlation lengths for other parameters. This gives an effective length of the chain of ~ 1 x 10^, 
which in turn gives statistical uncertainties in the parameter errors of about 1/v 10"^ = 1%, i.e., 
negligible. These uncertainties in the errors are consistent with the results from running multiple 
test chains. 

By making steps of the same size in the positive and negative directions for these parameters 
in constructing the Markov chain, we would implicitly assume that our prior knowledge of these 
parameters is a uniform distribution. This is not an accurate representation of our prior knowl- 
edge since, for example, binaries with periods between 10^ and 10^ years are not 10 times more 
common than binaries with periods between 10 and 10^ years. We therefore applied a prior to the 
likelihood function in the MCMC fitting where P and a are distributed evenly in logarithm and 
that the parameters ecos(a;) and esin(ti;) are uniformly distributed, rather than e and uj, to save 
the algorithm from unnaturally preferring circular solutions. This is equivalent to the /(e) = 2e 



^''For visual binaries, there is a 180° ambiguity in determining Q,, which can only be resolved with radial velocity 
information. 
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distribution as discussed by, e.g., iDuquennov &: Mavoij (| 19911 ). The very small effect of the choice 
of prior is discussed below P^l 

As a consistency check, we also ran our M CMC fitting code on the astronietric data for the 



binary L dwarf 2MASSW J0746+2000AB from iBouy et al.l (|2004l ). We found excellent agreement 



between the orbital parameters derived by us (using MCMC) and by Bouy et al. (using a variety 
of chi-square minimization approaches). Not only do the results agree to within the quoted errors, 
there is better than ~1% agreement on the best-fit results and better than ~20% agreement on 
the 95% confidence intervals. 



3.3.2. Fitting Results 

Figures [5] show the resulting probability distributions for the orbital parameters from the 
MCMC chain. The probability distributions are clearly not gaussian. For a given parameter, we 
adopt the median as the result and describe a confidence limit of X% as simply the ^ ± ggg bounds of 
the sorted sample. At 68(95)% confidence, we find a modest eccentricity of 0.25_g'-^gL'2Q|, an orbital 

period of lh.\_^^,^A yr and a semi-major axis of l'''l_-|^3L7) ™^s {'^■'^_Q2(Qi) ^U including the 
uncertainties in the plate scale and parallax). Two of the orbital angles are very well-constrained, 
the inclination i = 84.3_QgQ'yN deg (nearly edge-on) and the PA of the ascending node O = 

13.0_p'2L'gx deg. The final results are summarized in Table HI 

Figure [6] shows the strong correlation between the determination of the orbital period and the 
eccentricity. It illustrates that there are two classes of possible orbits: one branch having shorter 
periods and smaller semi-major axes and the other branch having longer periods and larger semi- 
major axes. Figure [7] shows that the short-P branch orbits have just passed apoastron [uo ~ 179°, 
~ 13°). This is the favored solution, with 98% of the steps in the MCMC chain residing in this 
branch (using e = as the dividing criteria in the P — e plane). However, a nearly circular orbit 



^*Soine care is also needed in handling the astrometric calibration errors of the three instruments used in our 
analysis. In practice, most of the constraints come from the fl'ST/WFPC2 discovery epoch from August 2000 and 
the six Keck LGS epochs from 2005-2008. Given the large time difference, the two datasets essentially constrain 
different portions of the parameter space. Some extra information is provided by the HST/ACS data but to a much 
lesser degree since it is contemporaneous with the Keck LGS data and has 2-4x larger errors. (In fact, the MCMC 
fitting gives basically the same results if the ACS data are excluded.) Thus simply applying the NIRC2 astrometric 
calibration error of 0.1% in pixel scale and 0.07° in PA to all six Keck epochs would incorrectly treat this as a 
random error, when in fact the uncertainty in the NIRC2 calibration globally impacts the overall solution of, e.g., 
the semi-major axis and the PA of the line of nodes. Therefore, we do the following; (1) we apply the errors in 
the NIRC2 calibration in quadrature to the HST/ACS and _FfST/WFPC2 measurements before orbit fitting, and (2) 
once the semi-major axis distribution has been determined, we apply the uncertainty in the NIRC2 pixel scale in 
determining the error on the total mass. The net effect is negligible, given the much larger ACS errors compared to 
the Keck LGS data and the fact that the error in the total mass is dominated by the parallax error (which is 15 x 
larger than the NIRC2 pixel scale uncertainty). 
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means it can be difficult to clearly distinguish whether the system has just passed apoastron or 
periastron, and thus a minority of the MCMC steps (2%) fall into the long-period branchP^I 

The MCMC fitting provides probability distributions for the orbital parameters, but does not 
provide a single best-fitting orbit per se, since a range of possible orbits fit the data with similar 
X^ values. One illustration of this is the result for lo = 179_-^Lgo-, deg, where the 95% confidence 
limits are broad enough to span both the short-period and long-period solution branches. Thus to 
plot orbits on the sky, we employ gradient-descent methods to find the best-fitting orbit with the 
MCMC-derived values as the starting point. Figures [9l and [TOl shows the resulting orbit, which has a 
period of 15.2 yr, a total mass of 0.0556 Mq and y^ = 0.9. To illustrate how the uncertainty in the 
orbital period impacts the orbit, we also show the best-fitting orbits found when fixing the period to 
12 and 20 yr, which have total masses of 0.0523 and 0.0590 Mq and x^ of 1.1 and 1.0, respectively. 
All three orbits reside in the short-period branch and show that the projected separation is now 
rapidly decreasing. The system is expected to be well-resolved again in the year 2011. 

Applying Kepler's Third Law to the period and semi-major axis distributions gives the poste- 
rior probability distribution for the total mass of the binary, with median of 0.0556 Mq, a standard 
deviation of 0.0018 Mq (3.2%), and a 68(95)% confidence range of about ±0.0018(0.0037) M© 
(Figure [TT|) . However, the MCMC probability distribution does not include the uncertainties in 
the parallax (1.6%) and the NIRC2 pixel scale (0.11%). By Kepler's Third Law, the quadrature 
sum of these errors amounts to an additional 4.9% uncertainty on the derived total mass. Since the 
MCMC-derived mass distribution is asymmetric, we account for this additional error in a Monte 
Carlo fashion; for each step in the chain, we draw a value for the pixel scale and parallax from a 
normal distribution and then compute the total mass. The resulting mass distribution is essen- 
tially gaussian (Figure [TT]) . Our final determination of the total mass is 0.056 it 0.003(0.006) Mq at 
68(95)% confidence. Thus, the total mass of this system is well-measured, with the parallax error 
being the dominant uncertainty. This is the coolest and lowest mass binary with a dynamical mass 
determination to date. 



3.3.3. Alternative Orbit Fits 

The WPFC2 discovery epoch in 2000 is obviously a key component to fitting the orbit. As 
described in § 12.2.11 we independent ly analyzed this dataset and greatly reduced the measurement 



errors compared those reported by iBurgasser et al.l (l2003bl ) . To examine the impact of this im- 



15 



Our inferred orbital period distribution is significantly longer than the 4- year estimate of lBurgasser et al.1 l|2003br ). 



Their original estimate was based on the proje cted physical separation at the discovery epoch, an assumed total mass 



of 0.07 Mq, and the statistical estimate from iFischer fc Marcvl (1199^) that the true semi-major axis is on average 



1.26x larger than the projected separation. The large discrepancy with our orbital period determination is not 
surprising, since a single epoch of imaging prov ides a highly uncertain period estimate — the likely true periods can 



surpnsmg, since a single epocn ot imaging prov iaes a nigniy 
span a factor of several greater or smaller (e.g. lTorreslll999l ) 
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provement, we also tried fitting the orbit using the original Burgasser et al. astrometry. Without 
the 2008 Keck data, our improved WFPC2 astrometry is essential for a well-constrained fit. How- 
ever, with the complete dataset, the fitted orbital parameters and the total mass are insensitive to 
the specific choice of WFPC2 astrometry, changing by less than la. 

Also, to check the effect of our assumed prior on the MCMC fitting (flat in logP and log a), 
we tried three alternative priors from the literature: 



1. Solar-type stars: iDuquennov &: Mavoij (jl99ll ) analyzed a well-defined sample of 164 nearby 
solar-type stars (spectal type F7 to G9) and found a log-normal distribution in orbital period: 



dN 



exp 



d{\ogP) 
where P is the period in days and diogp = 2.3 



-(logP-4.^ 
2<P 



(3) 



2. Ultracool visual binaries: lAllenI (J2007l ) conducted a detailed analysis of published imaging 
surveys of 361 ultracool field objects to model the separation distribution as a log-normal 
distribution: 



dN 



d{\og a) 



exp 



-(loga-0.86)^ 



(4) 



where a is the semi-major axis in AU and crioga = 0.28cj 



3. Ultracool visual and spectroscopic binaries: iMaxted &: Jeffries! ()2005l ) analyzed a sample of 
47 ultracool binaries with multiple radial velocity measurements and adopted a log-normal 
distribution truncated at large separations to match the known visual ultracool binaries: 



dN 



exp 



-(loga- 0.6)2 



2a2 

log a 



for a < 10 AU 



d(log a) 
where a is the semi-major axis in AU, cioga = 1.0, and dN/d{\oga) = for a > 10 AU. 



(5) 



^^Since this distribution is derived from analysis of ultracool field dwarfs, one might arguably choose this as 
the default prior. However, Allen's input data are restricted to imaging and do not include any spectroscopic 
binaries; therefore, the smallest semi-major axes are poorly constrained. For instance, the Allen distribution predicts 
basicall y no binaries at 1 AU or smaller separations ( > Scr events), perhaps at variance with spectroscopic binary 
studies dlVIaxted fc Jeffrieiboosi : iBasri fc Reinerslbood '). The young eclipsing M6.5-I-M6.5 binary 2MASS J0535-05 
( Stassun et al.ll2006l ) is also highly anomalous with this distrib ution (^ 8(7 even t ). Fina lly, the Allen analysis does not 
separately consider objects of different spectral types, whereas iBurgasser et al.l (|2006cr i have suggested the separation 
distribution for T dwarfs may be tighter than for the L dwarfs. Thus our default prior of flat in log a is a conservative 
choice. Note that longer period of 2MASS J1534— 2952AB we find relative to the original estimate of ~4 years 
(a ~ 1.1 AU) is in accord with the distribution proposed by Allen, which would indicate that such a ~1.1 AU system 
would be uncommonly rare (~ Sct) relative to the ~100 known ultracool binaries. 
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Figure [12] shows the posterior probabihty distributions for the total mass from the differeii t 
priors. Overall, the choice of prior has very little effect on the mass determination. The lAllenI (|2007l ) 
distribution favors slightly higher masses, but its results are consistent with the other priors. 

As a final indepe ndent check, we also fit our astrometry using the linearized least-squares 
fitting routine ORBIT (JForveille et al.lll999l ). using the MCMC-derived parameters as the starting 
guess. The ORBIT results are given in Table [H with a resulting x^ = 0.9, and show excellent 
agreement with the MCMC results. 



4. Discussion 



A primary goal of measuring fundamental properties for ultracool binaries is to compare the 
measurements against theoretical models of their physical properties. A number of studies have been 
published for the previous ultracool visual binaries with dynamical masses (§ 1), with subtle and/or 
overt differences in the ways that observations are compared to models. In the analysis that follows, 
we strive to clearly elucidate the comparison of our 2MASS J1534— 2952AB observations to the 
models, both in terms of its approach and limitations. From the standpoint of the observations, we 
have high quality measurements of (1) the total mass of the system and (2) the individual absolute 
magnitudes, with (3) the individual bolometric luminosities only slightly less reliable. (We have 
not measured the complete spectral energy distribution but have accounted for this uncertainty in 
computing Lboi in § 13. 2p . We now examine what can be learned from these data in concert with 
evolutionary models and theoretical atmospheres. 



4.1. Substellarity 

The most immediate result from our measurement is that 2MASS J1534— 2952AB is a bona 
fide brown dwarf binary. The total mass of 0.0 56 ifc 0.003 Mp) is below the solar-metallicity stel- 
lar/s ubstellar boundary of wO. 070-0. 074 Mq (e.g. lHavashi fc Nakandll963l : lKumarlll963l:lBurrows et al 
200l|), with the boundary increasing to higher masses for lower metallicities (jSaumon et al.lll994l ) 
— therefore, the individual components are clearly substellar. This is the second binary where 
both components are direc tly confirmed to be brown dwarfs, after the young eclipsing M6.5+M6.5 
binary 2MASS J0535— 05 ([Stassun et al.ll2006l ). and is the first such field binary in this category. 



4.2. Age 

Brown dwarfs follow a mass-luminosity-age relation. We have measured two of these quantities, 
the (total) mass and the luminosity, and by using evolutionary tracks we can deterinine th e third 
quantity, the age of the system. We use models from the Tucson group (JBurrows et al.lll997l ). which 



21 



provide predictions for Lboh and the "COND" models from the Lyon group (JBaraffe et al.ll2003l ). 
which predict both Lboi and absolute magnitudes. We conservatively assume the system is coeval 
and that the system is a true binary, not a partially (un)resolved higher order multiple system 

(§EID. 



For each tabulated model age, we use the individual absolute magnitudes and/or bolometric 
luminosities to calculate the mass of the components and then sum the masses. We then apply 
the observed total mass to determine the age range of the system. All measurement errors in 
Lboi and the total mass are accounted for in a Monte Carlo fashion, namely we repeat the model 
calculations over multiple realizations for the Lboi and the total mass values. We take great care 
to account for the covariance between the relevant quantities in the calculation. For instance, the 
total mass of the system and the luminosity both depend on the parallax, and thus their errors 
are positively correlated; we therefore draw the parallax values from a normal distribution and 
incorporate these in determining the Monte Carlo distribution of total masses and luminosities, 
which themselveed are then propagated in the model-based calculations. This approach results in 
a probability distribution for the system's age (as well as the other resulting parameters discussed 
below), which we summarize with the median value and confidence limits. 

Figure [13] shows the results of these calculations to determine the age of the system. For 
a consistent comparison between the Lyon and Tucson models, we use only the results derived 
from the Lboi measurements. However, the Figure also shows that using the absolute magnitudes 
predicted by the Lyon models would give similar results. 

We determine an age of 0.73ib0.07(0.15) Gyr from the Burrows models and 0.83ib0.08(0.18) Gyr 
from the Baraffe models at 68(95)% confidence. To construct a representatve "average" of the model 
results, we merge the results of the individual Monte Carlo calculations into a single distribution 
and compute its confidence limits. Thus, we assign an age of 0.78 it 0.09(0.18) Gyr (Table[5]). This 
is relatively youthful compared to the main-sequence s tars in the solar neighb orhood, e.g., ~95% of 
nearby solar-type stars have age estimates of >1 Gyr (JNordstrom et al.ll2004l ). However, the mean 
age of T dwarfs is expected to be younger than for field stars, since the known census is magnitude- 
limited and younger objects are b righter. The age d istribution of field ultracool dwarfs has been 
modeled by iBurgasserj (J2004l ) and lAllen et al.l (J2005l ) ; they generally find that field T dwarfs can 
span younger ages than for low- mass stars, though the predicted age distributions for both types 
of objects have large spreads. 

Kinematics provide an independent (albeit indirect) indicator of age, as older objects are 
expected to generall y show larger space motions due to their accumulated history of dynami- 
cal interactions (e.g. I WielenI EjjTtI ) . The tangential velocity of 2MASS J1534-2952AB {Vtan = 
17.3 lb 0.4 km s~^) is the second smallest measu red for T dwarfs, with only the T5.5 dwarf 



2MASS 1546-3325 being smaller (12.1±0.4 km s~^: lTinnev et al.ll2003l '). This is generally in accord 
with the 0.78±0.09 Gyr inferred from the evolutionary models, namely that 2MASS J1534-2952AB 
is among the youngest members of the nearby field population. However, since the measured Vtan 
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distribution of field T dwarfs is quite broad, with an unweighted aver age of 38.4 km s ^ a nd a 



standard d e viatio n of 20.4 km s ^ among the 21 unique objects in the iTinnev et alj (J2003l ) and 



Vrba et al.l (|2004l ) parallax samples, 2MASS J1534— 2952AB does not appear to be anomalously 
young for a field object. A radial velocity measurement is needed to determine the binary's space 
motion and thus better constrain its kinematics. 



4.3. Temperatures and Surface Gravities 

With the age of the system determined above, the combination of the observations and the 
evolutionary models provide highly precise values for the remaining physical parameters. The 
results derived from the two sets of evolutionary tracks are given in Table [5] and are computed 
from the same Monte Carlo approach that accounts for the covariance in the measurements. The 
Tucson and Lyon models give consistent values, with the Tucson models giving slightly larger 
radii and thus slightly cooler temperatures. Again, to compute a representative "average" for 
each parameter, we merge the Monte Carlo distributions computed from each set of models and 
compute confidence limits for the aggregate. This is not intended to be physically meaningful, 
but rather is a quantitative representation of the results that accounts for non-gaussian and/or 
inconsistent distributions from the two sets of models. We thus find radii of 0.0986 it 0.0015 Rq 
and 0.0993 ± 0.0017 Rq, effective temperatures of 1028 ± 17 K and 978 ± 17 K, and surface gravities 
of log{g) = 4.91 lb 0.04 and 4.87 it 0.04 for components A and B, respectivelvr^l 

To re-iterate, these properties are derived using only the measured total mass and resolved 
magnitudes/luminosities, along with the assumption that the system is coeval and composed of two 
components. No additional assumptions have been made to determine the individual masses. We 
also have avoided using spectral types and/or effective temperatures in this aspect of our analysis 
(contrary to some previous studies), as these quantities can introduce additional systematic errors 
and/or circular reasoning. For instance, it would be in correct to e n iploy the relations between 



spectral type and T^g from iGolimowski et al.l (|2004al ) or IVrba et al.l (J2004l ) to determine T^s for 



the two components and then compare to evolutionary tracks, as the Vrba et al. and Golimowski 
et al. relations are derived from the radii of field brown dwarfs predicted by the evolutionary tracks 
themselves. Likewise, it is not necessary to use Tcfj determinations from atmospheric models to 
determine the age of the system or the physical properties of the individual components in our 
approach. (See also § 14.51 ) 

As already noted above, our Monte Carlo calculations account for the covariance in the mea- 
surements. One important effect is that the uncertainties derived from the resolved magnitudes and 
luminosities of the two components are correlated, since they all depend on the measurement uncer- 



^^Since all these derived properties rely on the model-derived age, their distributions from the Monte Carlo calcu- 
lations are fairly correlated. However, the formal confidence limits on the quantities are so small that this correlation 
is unlikely to be significant for any future analyses. 
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tainties in the integrated-light photometry of the system. As a consequence, the relative temperture 
difference between the two components (ATefj = SOJ^iq K) can be calculated to higher precision 
than would be indicated by the uncertainties in the individual T^s det erminations (v2 x 17 K = 
24 K). This ATgg agrees with the ~70 K difference expected from the iGolimowski et al.1 (J2004al ) 
polynomial fits for the 0.5 subclass difference between the two components. (The ATeff from the 
Golimowski et al. fits would be about twice as large for a 1 subclass difference.) 

The physical parameters for 2MASS J1534— 2952A and B are in general agreement with pre- 
vious determinations for the properties of field T dwarfs. However, our values have much higher 
precision, because the accurate total mass measurement leads to a small age range, which leads to 
strong constraints on the radii and thus small uncertainties on the derived T^n and log{g) values. 
(We discuss this further in § 14.51 ) T dwarf surface gravities have been inferred to be l og(q) = 4.5—5.5 
by comparing theo retical model atmos pheres to optical spectra (JBurrows et al.ll2002), near- IR colors 



2nn6a 



Saumon et al. 



2006 



and line strengths (IKnaPD et al.ll2004l). arid low-resolution near-IR/in id-IR spectra (JBurgasser et al. 



Leggett et al. 



2007a 



Gushing et al.l 120071 ). This range encompasses our 



log(g) determinations for 2MASS J1534-2952AB. 

On the other hand, our very precise temperatures for the two components of 2MASS J1534— 2952AB 
are discrepant with previous studies of T dwarfs: the Tefj's of 2MASS J1534— 2952AB appear to 
be cooler than determined previously for mid-T dwarfs. These discrepancies occur for two separate 
comparisons. 



4-3.1. Temperature discrepancy with evolutionary models 

Temperatures for field T dwarfs have been inferred by combining accurate Lboi determinations 
with radius predictions from evolutionary tracks. This approach is expected to be reasonably 
accurate, since th e radii of brown dwarfs older than ~100 Myr are predicted to vary by <30%. 
Golimowski et al.l (|2004al ) ad opt a typical age of 3 Gyr and a plausible range of 0.1-10 Gyr in 
cor nputing Tpg frorn Lboi, and lVrba et al.l (J2004l ) adopt a radius range of 0.90 ±0.15 i?jup predicted 
by iBurgasserJ (120021 ) simulations of the solar neighborhood assuming a constant star formation 
history. (Both studies use the Tucson evolutionary tracks.) There are 8 T4.5-T6.5 dwarfs in these 



studies (after updating the spectral types to the latest classification by iBurgasser et al.l (l2006bl'l) 



and all of th ese appear to be single based on high angular resolution imaging (IBurgasser et al 



2003bl . l2006d : Liu et al., in prep). We compute the average and standard deviation of both studies 
to obtain T^s = 1216 ± 20 K for T4.5 (2 objects), 1146 K for T5.0 (interpolated), 1077 K for T5.5 
(1 object), 1014 ±33 K for T6.0 (3 objects), and 950 ±106 K for T6.5 (3 objects). We have excluded 
the T6±T8 binary 2MASS J1225-2739AB and the pecuhar T6 dwarf 2MASS J0937±2931, and we 
have assumed the T4.5 dwarf 2MASS J0559— 1404 is an equal-magnitu de unresolved binary based 
on its pronounced overluminosity (e.g., Figure 3 of iBurgasserl l2007al ). Neither sample contains 
any T5.0 objects, but the 4 other subclasses almost exactly follow a straight line, so we linearly 
interpolate to find Tcs for T5.0. 
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In comparison to the field objects, the components of 2MASS J1534— 2952AB appear to have 
~100 K cooler temperatures relative to their spectral subclass (~ 120 it 35 K for the primary and 
~ 100 lb 35 K for the secondary, where we have adopted a 30 K uncertainty for the T5 and T5.5 field 
objects based on the other subclasses with more than one object). The disagreement is modest, 
and a more definitive comparison is hampered by the few Tf.Q determinations (i.e., parallaxes) for 
T4.5-T6.5 dwarfs. Nevertheless, the result is potentially intriguing. 



In particular, iMetchev &: Hillenbrand! (|2006l ) have noted perhaps a similar effect for the three 
known late-L (L7-L8) dwarf companions to field stars. More precise T^s estimates can be obtained 
from evolutionary models for these companions than for field objects by incorporating the age 
estimates of their primary stars. (See § 14.51 for details.) Metchev & Hillenbrand find that the 
L dwarf companions appear to be ~100-200 K cooler than single field late-L dwarfs. They raise 
the possibility that the model-derived radii are at fault, either due to incorrect cooling rates or 
systematic overestimate of the field dwarf ages. However, they prefer the hypothesis that the 
discrepancy is a manifestation of an unanticipated surface gravity dependence of the L/T transition, 
causing younger L/T transition objects to have cooler temperatures than older ones. This is 
motivated by their analysis of the young (0.1-0.4 Gyr) L7.5 companion HD 203030B and apparently 
supported by the young (0.1-0.5 Gyr) T2.5 compan ion HN Peg B, which also appears to be R i200 K 



cooler than field objects of the same spectral type (JLuhman et al. 



2007 



Leggett et al.ll2008l l. 



We find that the T5.0 and T5.5 components of 2MASS J1534-2952AB may also be «100 K 
cooler than comparable field objects. Since these two components are later-type than the L/T 
transition (e.g., their positions in IR color- magnitude diagrams is coincident with the locus of 
mid/late-T dwarfs with blue near-IR colors), this may suggest that the Tgfj discrepancy might not be 
solely associated with the L/T transition. Instead, 2MASS J1534— 2952AB and the aforementioned 
L/T companions may indicate that a systematic error in the estimated ages and radii of field late-L 
and T dwarfs is the culprit. 

The sslO% temperature discrepancy for 2MASS J1534-2952AB amounts a ?a20% underes- 
timate of the radii. For a fixed value of Lboi (which is the appropriate constraint here) , the 
Burrows et al.l (120011 ) scaling relations give 

8.56 



tr^ K 



(6) 



where t is the age and R is the radius. This agrees well with the exponent value of 8.2-8.4 
extracted from the Tucson models for sources of log(Lboi/-^o) = —5.0 to —5.1. Propagating the 
±3% uncertainties in the T^g disagreement, the implied age overestimate is a factor of 6 it 3, meaning 
implied ages of ~0. 3-1.0 Gyr for the field populationPi 

The same discrepancy can be seen in an alternate fashion, namely by comparing the lumi- 
nosities for the same T4.5-T6.5 field objects: log(Lboi/-^o) = -4.79 ± 0.04 for T4.5, -4.89 for 



^** Applying the same scaling relation to the 10-15% radius discrepancy found bv lMetchev fc Hillenbrand! (200a) for 
objects at the L/T transition implies a factor of 2-3 overestimate in the representative age of the field population. 
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T5.0 (interpolated), -4.99 ± 0.01 for T5.5, -5.17 ± 0.29 for T6.0, and -5.22 ± 0.20 for T6.5. 
The luminosities of 2MASS J1534— 2952AB are comparable to the field objects of similar type 
(—5.015 lb 0.019 for T5.0 component A and —5.093 it 0.019 for T5.5 component B). Thus, in order 
for all the objects to have similar temperatures and Lboh they must have about the same radius 
and thus about the same age as 2MASS J1534— 2952AB. In other words, the measured total mass 
of 2MASS J1534— 2952AB is too small (by a factor of ~2) compared to the mass expected from 
the evolutionary models for 3 Gyr objects with log(Lboi/-^0) ~ —5.0. 

A representative age of ~0.5 Gyr for the field population is not ruled out given the state of 
the observations. Though Golimowksi et al. did consider the range of 0.1-10 Gyr, they adopted 
a nominal age of 3 Gyr in determining T^s for field dwarfs, based on th e 2-4 Gyr age estimate 
from the tangential velocities of ultracool dwarfs by iDahn et al.l ()2002l ). A younger age could 
be accomodated, since the tangential velocity of a population is only an approximate statistical 
estimate of its age. Indeed, kinematic analysis of the space motions of L and T dwarfs suggests a 
younger age of ssO.5-2 Gyr (lOsorio et al.ll2007l ). Similarly, the radii of 0.90 ±0.15 i?jup adopted by 
Vrba et al. is based on a mass function where dN/d M ^ M~^\ a soin ewhat st eeper mass function 



would lead to younger typical ages (e.g.. Figure 8 of iBurgasserl |2004| . though iMetchev et al.l 12007 
suggest dN/dM ~ M^ based on a small sample of T dwarfs). Thus, the discrepancy of evolutionary 
model-derived temperatures between objects of known mass/age and the field population can be 
plausibly explained by a modest overestimate of the ages of the field population. A larger sample 
of ultracool dwarfs with known masses and/or ages is needed to better explore this issue (§ 14. 5p . 



4-3.2. Temperature discrepancy with model atmospheres 

The spectrum of 2MASS J1534— 2952AB has not yet been fitted with model atmospheres due 
to its composite n a.ture. In fact, the T dwa rf class as a whole has not been extensively subjected to 
such comparisons. iBurgasser et al.l (|2006al ) determined Tes for a sample of sixteen T5.5-T8 dwarfs 
by comparing near-IR spectral indices to condensate(dust)-free atmosphere models from the Tucson 
group. They determined T^s = 1020 — 1100 K for one T5.5 dwarf, and an unweighted linear fit 
of atmosphere-derived T^g versus spectral type for their sample (excluding the peculiar T6 dwarf 
2MASS J0937+2931) gives 



TcS, atmosphere = 1090 - 126 X {SpT - 25.5) 



(7) 



where SpT = 25.5 for T5.5, SpT = 26 for T6, etc. The RMS about the linear fit is 50 K, which 
we adopt as the uncertainty (a value somewhat larger than the ±10 K to ±40 K range computed 
for individual objects in their sample). Extrapolating the linear fit gives Tefj = 1160 K for T5. 
This is obviously approximate, e.g., given the potential systematic effects in the models and the 
spectral classification scheme, though this value agrees with the T^q = 1150 — 1200 K found by by 
fitting model atmosp heres to the 0.95-14.5 fxm spectrum of the T4.5 dwarf 2MASS J0559— 1404 
(JGushing et al.ll2007l ). 
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Thus, model atmospheres indicate T^s = 1160±50 K and 1090±50 K for 2MASS J1534-2952A 
and B, respectively. The temperatures we find using evolutionary tracks appear to be cooler by 
~100 K at modest significance (—130 ± 50 K for component A and —110 it 50 K for component B, 
if we assume that the errors add in quadrature). We cannot objectively discern if the problem 
lies in the evolutionary tracks, the model atmospheres, or bot h. However, the evolu tionary models 
are thought to be robust to the principal input uncertainties (IChabrier et al.ll2000l ). On the other 
hand, the model atmospheres are quite uncertain. Even though the spectral appearance of mid- 
and late-T dwarfs is relatively simple — dominated by collision-induced H2, H2O, and CH4 in 
the near-IR and the wings of the K I 0.77 fira resonance line in the far-red — the line lists for 
H2O and CH4 are known to be incomplete, and the input physics to the atmosphere models are 
complex. Current a tmospheres generally match the observed spectra of late-T (T6-T8) dwarfs, but 



not ex actly so (e.g. iBurrows et al. 



2006 



Burgasser et al.ll2006al : ISaumon et al. 



2006 



Leggett et al. 



2007al ) r^l Therefore, while a larger sample of objects is needed both for dynamical mass determi- 
nations and model atmosphere fitting, the plausible hypothesis is that the observed discrepancy 
arises from an overprediction of T^s by current model atmospheres. 



4.4. Color-Magnitude and Hertzsprung-Russell Diagrams 



We have directly measured the total mass of the 2MASS J1534— 2952AB system. However, 
using the evolutionary tracks to determine the physical properties also implicitly determines the 
mass ratio, since the model-derived age and observed individual luminosities translate into individ- 
ual masses (again with the assumption that the system is composed of only two components). We 
infer the mass ratio of the system from the ratio of the bolometric luminosities, as this is very ro- 
bust. TojlhistratethiSj_consider the analytic scaling relation for solar-metallicity substellar objects 



from Burrows et al. 



(2001 



Lb. 



ol 



^2.64 ^-1.3 ^a35 



(8) 



where M is the mass, t is the age and kr is the Rosseland mean opacity. We measure a 0.078 it 
0.016 dex difference in Lboi between the two componentCj which leads to a mass ratio q = 
Mb/Ma = 0.934 it 0.007. The uncertainty in the mass ratio is small due to the weak depen- 
dence of mass on luminosity at fixed age. Using the actual tabulated Tucson and Lyon models 
and again keeping careful track of the covariance in the calculations, we compute a final value 
of (/ = 0.936lo'oo8' where the error includes the uncertainties in the model-inferred age and the 



'^^For instance, recognizing these limitations, iBurgasser et al.l (2006a|) chose to cahbrate the model atmosphere 
predictions empirically using the well-studied T7.5 dwarf Gl 570D in fitting models to late-T dwarf spectra, as 
opposed to using the atmospheres directly. 

^"At face value, the results in Table |3] would give an uncertainty of \/2 x 0.018 — 0.025 dex in the Lboi difference, 
but this would include the uncertainties in the distance modulus (0.04 mag) and the integrated-light photometry 
(0.03 mag), which are common to both components. 
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observed Lboi difference. (The Tucson and Lyon models give basically identical results for q.) 
This gives individual masses of 0.0287 ± 0.0016 Mq (30.1 ± 1.7 Mjup) and 0.0269 ± 0.0016 Mq 
(28.2 lb 1.7 Mjup) for components A and B, respectively. 

We first compare the individual components against the COND evolutionary models of the 
Lyon group, which provide predictions for the absolute magnitudes and colors. The model predic- 
tions are generated for the CIT photometric system, so we tr ansform our resolved MKO photometry 
for 2MASS J1534-2952AB to this system using the results of lstephens fc LeggettI ^M)- FigureH 
shows that the models are somewhat too red compared to the data. This is not surprising, as model 
atmospheres for T dwarfs are known to be deficient in the CH4 and H2O opacities relevant at these 
wavelengths (e.g. iLeggett et al.ll2007al ). The plotted COND models are also computed only for 
solar-metallicity, and a non-solar metallicity for 2MASS J1534— 2952AB would impact the colors 



and magnitudes (e.g. iLiu et al. 



2006a 



Burgassen l2007bl ) . Indeed, current models do not exa ctly 



match the observed color-magni tude loci for field T dwarfs (e.g.. Figure 8 of lKnapp et al.ll2004l and 
Figure 7 of Burrows et ahlbood ). Nevertheless, 2MASS J1534-2952AB will provide a strong test 
to for future models, since the components' magnitudes, colors, and masses are very well- measured. 

With the individual mass estimates and an independent determination of T^a, it is possible to 
directly test different evolutionary tracks using the Hertzsprung-Russell (H-R) diagram. We use the 
values of T^s = 1160 ±50 K and 1090 ±50 K for 2MASS J1534-2952A and B, derived in §[43]from 
model atmosphere studies. Figure [15] shows the individual components on the H-R diagram and 
compares these to the Tucson and Lyon evolutionary tracks. The locations of the two components 
disagree with both sets of models (which agree very well between themselves) . 

The likely interpretation is that the temperatures from the model atmospheres place the com- 
ponents to the left of the evolutionary tracks, i.e., too warm. As discussed in § 14.31 the model 
atmospheres are a significant source of the uncertainty in placing the components on Figure [151 A 
possible systematic error of only ~100 K would be sufficient to resolve the discrepancy with the 
data. Therefore, while acquiring resolved spectra of the two components could help refine the tem- 
perature determination, the systematic uncertainties in the atmosphere models will still hamper 
accurate placement on the H-R diagram. We discuss this further in the next section. 

The opposite interpretation is that the evolutionary models are incorrect, leading to a ~50% 
overprediction of the luminosities and a ~20% overprediction of the radii, given the component 
masses. Equivalently, if one were simply to assume the H-R diagram positions are accurate, the 
evolutionary models would suggest individual masses of around 0.05 Mq and 0.06 Mq from the 
Lyon and Tucson models, respectively, i.e., nearly a factor of two overestimate in the masses. 
While it may be that the evolutionary models are so substantially incorrect, such a conclusion is 
not compelling at this point, given the plausible errors in the Tefj determinations!^^! 



^^Note that a systematic error in our model-derived mass ratio cannot resolve the H-R diagram discrepancy. Since 
the ratio is basically unity, correcting any errors in q would move one evolutionary track closer to one component, 
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Direct mass determinations for the individual components from radial velocity monitoring 
and/or absolute orbital astrometry will help to further characterize the system. Such data will 
directly test the q = 0.936^g'QQg determined from the evolutionary tracks. The expected maximum 
radial velocity difference of the two components is only 4.6 km s~^. Since the two components 
are nearly equal mass and brightness, the orbital motion will be very difficult to detect in the 
integrated-light spectrum. Resolved AO spectroscopy will be required, and the small amplitude 
will make it a challenging measurement giv en the few kin s~^ accuracy that has been achieved for 
T dwarfs on the largest existing telescopes (jOsorio et al.l 120071 ). 



Individual mass measurements can in principle also test the evolutionary tracks directly. One 
can estimate the age of each component from its mass and luminosity (as we have done using 
the total mass) and see if the ages indicate coevality for the system. However, given the near- 
equal flux ratio of this system (and most ultracool binaries), this coevality test is unlikely to be 
very discriminating. Moreover, individual masses cannot resolve the discrepancy seen in the H-R 
diagram (Figure [T5]) . which largely arises from the uncertainties in the model atmospheres. 



4.5. Future Tests of Theory with Field Substellar Binaries 

With the advent of LGS AO on the largest ground-based telescopes, we can expect an in- 
creasing number of dynamical masses for ultracool field dwarfs in the near-future. The most 
useful systems for testing theory will be those with both independent mass and age determina- 
tions, namely binaries that are associated with open clusters/groups and/or field stars of known 
age. The former w ill present a significant technical challenge, e.g., ultracool binaries in the Hyades 
{d = 46.3 ± 0.3 pc; IPerrvman eiTaPllQQsl ) and Pleiades {d = 133.5 ± 1.2 pc; ISoderblom et ahlbood ) 
with suitably short orbital periods are unresolvable with current technology and thus none are 
currently known. Ultracool binary companions to field stars are extremely rare and thus while 
very valuable systems, these will only probe a very limited range of spectral type, age, and mass: 
only four systen is are known with suitabl y short orbital periods {P < 50 yr) — the T1-I-T6 b i- 



nary e Ind Bab (JMcCauehrean et al.ll2004l^. t he L4 +L4 binary H P 130948BC (|Potter et al.1 I2OO2I ') . 
the L4.5-FL6 binary GJ 417BC (IBouv et al.l l2003l: fcizis et al.l [20oil. and the L4.5+L4.5 binary 
GJ 100 IBC (JGolimowski et al.ll2004bl ). Therefore, there is significant motivation to develop analy- 
ses that employ masses derived from the much more numerous field binaries. In this regard and as 
illustrated by our analysis for 2MASS J1534— 2952AB, one can identify two orthogonal pathways to 
confront theory: (1) comparison to evolutionary tracks and (2) comparison to atmospheric models. 



while the other track would move farther away from the other component. The evolutionary models could be brought 
into agreement with the data for component A for q < 0.4, but this is implausible given the nearly equal magnitudes 
of the two components. Such a small q would also exacerbate the disagreement between the observations and the 
models for component B. 
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4-5.1. Comparison to evolutionary tracks ("H-R Diagram Test") 



Direct measurements of Lboh ^eflfi and mass (or age) for brown dwarfs enable use of the H- 
R diagram, by comparing the observations to evolutionary tracks that correspond to the measured 
masses of the objects. As illustrated by Figure [T5l the Lyon and Tucson tracks differ at the 5-10% 
level in mass, and thus mass determinations of 2-3% accuracy could discriminate between the two 
models, if Lboi and T^s can be well- measured. (Improvements in the parallaxes of many ultracool 
binaries will also be needed to achieve such accurate masses.) Accurate measurements for Lboi are 
largely straight-forward, as good as a few percent (e.g., § l3.2p . However, direct Tes determinations 
are extremely challenging, since radius measurements are needed. Brown dwarfs are too small and 
faint to be resolved with current or planned interferometers, and no eclipsing ultracool field binaries 
are yet known. Thus, T^s must be derived from modeling the observed colors, magnitudes, and/or 
spectra; the approach currently suffers from uncertainties at the level of a few to several hundred 
Kelvin and systematic errors that are difficult to quantify (e.g. lCushing et al.ll2007l ). In comparison. 
Figure [15] shows that T^g determinations good to <30 K are needed. Therefore, decisive tests of 
evolutionary tracks using field binaries will be challenged by this uncertainty in Teg , in the absence 
of direct radius measurements. 



4.5.2. Comparison to atmospheric models ("Age/Mass Benchmark Test") 

Brown dwarfs obey a mass-luminosity-age (M, Lboi,i) relation, and for most field objects nei- 
ther the mass nor the age is known. A commonly used approach to circumvent this limitation is to 
study brown dwarfs that are compani ons to main-sequence stars, where (indirect) age estimates are 
available from the primary star (e. g. 



2004 



Metchev &: Hillenbrand 



20061 : 



Saumon et al. 



Liu et al. 



200c 



2006a 



Burgassei 



2001 



McCaughrean et al. 



Geballe et al. 

2007bl ) . This approach can also be 



applied t o members of coeval clusters /groups and corn panions to post-main-sequence stars of known 



age (e.g. iKirkpatrick et al.lll999l : iPinfield et al.l 120061 ). In these situations, Lboi and t are known, 
and combined with evolutionary models, one can derive M and consequently T^s and log(5f). Then 
the observed colors, magnitudes, and spectra can test the accuracy of atmospheric models with 
the same Tefj and log(5f). Examination of the known "age benchmark" T dwarfs in this fashion 
finds that the properties deduced from atmospheric models are in good agreement with those from 
the evolutionary models, within the uiicertain ties in the ages and metallicities of the primary stars 
(|Burgasserll2007bl : iLeeeett et allboOTal . bood ). 



We suggest that, in an analogous fashion, field binaries with known masses can also serve 
as "benchmark" objects. In this case, M and Lboi are known, and combined with evolutionary 
models, one can derive t, as we have done in § 4.1. This provides T^g and log (5) and thereby 
allows tests of atmospheric models. The chain of analysis is identical to brown dwarf companions 
of known age: given independent knowledge of two quantities out of {M, Lboi, t}? use evolutionary 
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models to derive the thirdo In both cases, these benchmarks can also serve as anchor points for 
direct empirical calibra tion of spectroscopic diagnostics of log{g) and T^s (e.g. iPinfield et al.ll2006l : 
Burgasser et al.ll2006al ). 



At face value, using objects that are age benchmarks or mass benchmarks is less fundamen- 
tal than direct tests of the evolutionary models using the H-R diagram. However, in practice the 
Benchmark Test is much more feasible to implement and subject to much smaller systematic errors. 
In the absence of direct radius measurements, the H-R Diagram Test is held hostage to the system- 
atic errors in the determination of T^s from atmospheric models. In contrast , the Benchmark Tes t 
relies on the evolutionary models, which are thought to be more robust (e.g. IChabrier et al.l 12000 ). 
In short, given the choice of relying on atmospheric models (H-R Diagram Test) or evolutionary 
models (Benchmark Test), the evolutionary models are likely to be preferred. 

To assess the relative utility of "age benchmarks" (brown dwarf companions to stars) compared 
to "mass benchmarks" (brown dwarfs with dynamical masses), we turn to Equation[8l For an object 
with a measured Lboi and ignoring the weak dependence on kr, given a measurement of M or t with 
accompanying uncertainty of 6M or 5t, the fractional error in the remaining quantity is related by: 



t M 



(9) 



Typi c al uncertainties in the ages of main-sequence fie ld stars are about 50-100% (e.g. lKirkpatrick et al 



2001 



Liu et al. 



2006a 



Metchev &: Hillenbrandll2006l ). and thus age-benchmark objects would have 



a 25-50% uncertainty in the mass inferred from evolutionary modelso In contrast, dynamical 
masses of w5-10% accuracy will be possible in the next few years for many ultracool dwarfs, as we 
have already achieved with 2MASS J1534— 2952AB, and hence mass-benchmark objects will have 
only a ~ 10-20% error in the age inferred from evolutionary models. So overall, age benchmarks 
can be expected to have {5t/t, 5M/M} w {50-100%, 25-50%} while mass benchmarks would have 
uncertainties of order {5t/t,5M/M} w {10-20%, 5-10%}. 



We can use the analytic fits to evolutionary models from iBurrows et all ( 200ll ) to gauge the 
relative accuracy on Tefj and log{g) derived from both types of benchmarks. Using standard error 
propagation and assuming uncorrelated errors, we find for mass benchmarks : 



sz 



eflf 



6log{g) 



6M\'' / (JLboi 

0.180 + 0.246 — — 

^ M J \ Lboi 

1-63^^ + 0.088—^ 



(10) 
(11) 



^■^The analogy between mass benchmarks and age benchmarks is an imperfect one, since brown dwarf companions 
also have metallicity determinations from their parent star, whereas field binaries do not. However, current studies 
of age benchmarks largely rely on evolutionary models computed for solar metallicity. 

^''Barnes (2007) find that ages for solar-type stars derived from gyrochronology can have errors of only 15-20%. 
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VioooKy yV 0.05 J \ o.io J ^ ' 



W5M/mV ( 5L),o\/lY 



flogte) = 0.082 dex y ^^^^j + l^lOS^n^J <''> 

where (JTefj , 5M, (5Lboi j and (51og((7) are the uncertainties in the temperture, mass, luminosity, and 
surface gravity, respectively. And then for age benchmarks, we find: 



If - MhH-tJ <"> 



.St\ ( (5Lbol^ 



b\oz{g) = \/(0-734yl +(0-517-^) (15) 




2 / rr /r \ 2 



VioooKy y\o.5j \ o.io J 



(16) 



^t/t\ , /^n i/i'^-^bolZ-^t ^ 



0.5 ) 



Slogig) = 0.37 dex ^/ ( I^ ) + ( 0.14 ""-;^^^°' ) . (17) 



Thus with representative values for the fractional errors in age (50%), luminosity(10%), and mass 
(5%), we see that T^s and log{g) are better constrained by a factor of ~5 using mass benchmarks 
than age benchmarks. Figure [TBI plots the derived analytic estimates for both types of benchmarks. 
These contour plots provide a convenient means to gauge the expected errors in T^n and log(g) 
determinations from benchmarks. The morphology of the contours also illustrates whether the 
observational errors in the age, mass, and/or Lboi dominate the uncertainties in T^s and log{g). 
For the specific case of 2MASS J1534— 2952AB, there is good agreement between the analytic 
estimates and the values derived directly from the actual evolutionary models (Table [5]). 



5. Conclusions 

We have determined the first dynamical mass for a binary T dwarf, the T5.0+T5.5 system 
2MASS J1534-2952AB, by combining six epochs of Keck LGS AO imaging from 2005-2008 with 
three epochs of HST imaging obtained in 2000 and 2006. Both datasets achieve milliarcsecond 
accuracy or better for the relative astrometry of the two components and are validated through 
extensive testing with images of simulated binaries. We employ a Markov Chain Monte Carlo 
analysis to determine the orbital parameters and their uncertainties. The time baseline of our 
complete dataset covers about half of the total period. We find that the orbital motion of the 
binary is viewed in an almost edge-on orientation and has a modest eccentricity. Our determination 
of a 15.ll^'g yr orbital period is significantly longer than the original 4- year estimate, as by chance 
the binary was at a very small projected separation when discovered in 2000. 
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The total mass of the system is 0.056 it 0.003 Mq (59 it 4 Mjup), including the uncertainty in 
the parallax. This is the second brown dwarf binary directly confirmed, the first among the field 
population. It is also the coolest and lowest mass binary with a dynamical mass determination to 
date. 

With very accurate measurements of the total mass and the bolometric luminosity (log(Lboi/-^0) 
—4.751 lb 0.011), we use the Tucson and Lyon evolutionary tracks to determine the remaining 
physical properties for the system. The two sets of models give largely consistent results, which 
highlights the difficult of distinguishing between them even with such precise observational data. 
We average the model results to represent the final determinations. We find a relatively youthful 
age for the system of 0.79ib0.09 Gyr (lo"), consistent with its low tangential velocity relative to 
other field T dwarfs. The remaining physical parameters of the individual components are then 
fully determined: radii of 0.0986 ± 0.0015 Rq and 0.0993 ± 0.0017 Rq, effective temperatures of 
1028 ± 17 K and 978 ± 17 K, surface gravities of \og{g) = 4.91 ± 0.04 and 4.87 ±0.04, and masses of 
0.0287 ±0.0016 Mq (30.1 ±1.7 Mjup) and 0.0269 ±0.0016 Mq (28.2 ±1.7 Mjup) for components A 
and B, respectively. We take care to account for the covariances inherent in the measurement 
uncertainties, by using a Monte Carlo approach to derive these physical quantities from the evolu- 
tionary models. Our approach also assumes that the system is coeval and composed of only two 
components. 

These precise determinations for 2MASS J1534— 2952A and B are in general accord with the 
Teff and \og{g) values found previously for field T dwarfs based on model atmospheres and with 
the ages of T dwarfs predicted by Monte Carlo simulations of the solar neighborhood. However, 
upon closer scrutiny, there are two potential discrepancies with past studies. Both suggest that the 
temperatures of field T dwarfs may be overestimated by ~100 K, though we stress that the two 
discrepancies must arise from independent effects. (1) The temperatures of 2MASS J1534— 2952A 
and B appear to be cool er than field objects of compar able spectral type. T^^^ re sembles discrepan- 



cies previously noted bv lMetchev fc Hillenbrand! (|2006l ) and lLuhman et al.l (J2007l ) for late-L/early-T 



dwarfs that are companions to young main-sequence stars. They have hypothesized that the effect 
is due to the gravity sensitivity of the L/T transition. The fact that this discrepancy also occurs 
for 2MASS J1534— 2952AB suggests instead that the problem may arise from a factor of ~ 6 ± 3 
overestimate in the adopted ages of field objects when determining their temperatures using evolu- 
tionary tracks. Ages of ~0. 3-1.0 Gyr are preferred based on this binary. (2) The temperatures of 
2MASS J1534— 2952A and B are slightly cooler than inferred for other mid-T dwarfs from model 
atmospheres. Detailed analysis of the system's integrated-light and resolved spectra with model 
atmospheres is needed to directly assess the T^fj and log(5f) of the two components and to refine 
the comparison with the values derived from evolutionary models. 

The positions of the two components on the H-R diagram are discrepant with theoretical 
evolutionary tracks corresponding to their individual masses. In fact, taken at face value, using 
the H-R diagram positions to infer masses from the evolutionary tracks would lead to masses 
of ~0. 05-0. 06 Mq, about a factor of two larger than the actual measured masses. While this 
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discrepancy could stem from large systematic errors in the luminosities (~50% errors) and/or radii 
(~20% errors) predicted by evolutionary models, the likely cause is that temperatures from model 
atmospheres are too warm by ~100 K for mid-T dwarfs. This highlights the need for continued 
improvements to the model atmospheres. 

Future monitoring of 2MASS J1534— 2952AB will help to refine its orbit and its dynamical 
mass. The orbital separation of the system is now rapidly decreasing and will not be readily 
resolvable again until around 2011. At the same time, an improved parallax for the system will 
be required — the uncertainty in the total mass from the orbit fitting is 3%, compared to the 
5% that arises from the uncertainty in the parallax. Radial velocity monitoring and/or absolute 
astrometry will directly determine the individual masses and test if the system is a higher order 
multiple. However, given the very similar fluxes of the two components (implying nearly equal 
mass) , individual mass measurements are unlikely to resolve the discordant H-R diagram position of 
the two components relative to evolutionary tracks. This problem is likely driven by the systematic 
uncertainties in current model atmospheres for T dwarfs. 

The fundamental characteristic of the field population is that it spans a range of (largely 
unknown) ages. However despite this uncertainty, field brown dwarf binaries can strongly test 
theoretical models, if analyzed appropriately. These systems will be especially valuable in light of 
the current paucity of eclipsing field ultracool binaries and resolvable, short-period ultracool binaries 
in open clusters/groups CJ Specifically, attempts to directly test different evolutionary tracks by 
placing ultracool objects on the H-R diagram (the "H-R Diagram Test") will be challenging, given 
the similarity between the tracks and the difficulty in independently determining Tefj with model 
atmospheres. Instead, atmosphere models can be confronted against log(g) and T^s values for 
ultracool objects as derived from the evolutionary models, which can be exceptionally precise (the 
"Benchmark Test"). This approach has previously been applied to single brown dwarfs that are 
companions to stars of known age ("age benchmarks"). We suggest that in an analogous fashion, 
field ultracool binaries with dynamical mass determinations ("mass benchmarks") can test the 
model atmospheres. In fact, given the plausible observational uncertainties, mass benchmarks are 
likely to provide stronger constraints (by a factor of ~5) on log{g) and T^s than age benchmarks, 
since dynamical masses can be determined far more accurately than ages for main-sequence field 
stars. With the widespread advent of LGS AO on the largest ground-based telescopes, we can look 
forward to a rapid increase in dynamical mass determinations for low-mass field binaries and thus 
substantial advances in our understanding of the properties and evolution of substellar objects. 

We gratefully acknowledge the Keck LGS AO team for their exceptional efforts in bringing the 
LGS AO system to fruition. It is a pleasure to thank Antonin Bouchez, David LeMignant, Marcos 
van Dam, Randy Campbell, Al Conrad, Jim Lyke, Hien Tran, Robert LaFon, Kenny Graves, Cindy 



^''For reference, about 1 out of a 1000 stars are eclipsing binaries, whereas there are only ~600 L and T dwarfs 
known. 
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Fig. 1.— Images of 2MASS J1534-2952AB from HST/WFPC2, Keck LGS, and HST/ACS, ar- 
ranged chronologically in each column. Each image is 1.24" (16.8 AU) on a side, with the orientation 
indicated by the compass roses. Note that the sky directions for the ACS images are not orthog- 
onal, due to optical distortion in this instrument. (We chose not to rotate the HST images to the 
orientation of the Keck image for this figure, but instead to display the actual images so as to 
preserve the quality of the data.) The greyscale images use a square-root stretch. The contours are 
drawn from 80, 40, 20 and 10% of the peak pixel. Due to the effect of atmospheric dispersion, the 
Keck LGS J and -H-band images are slightly elongated in the vertical direction (which coincided 
with the elevation axis for these observations). 
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Fig. 2.— Top panel: Near-IR spectrum of 2MASS J1534-2952AB obtained with the IRTF/SpeX 
spectrograph. Other panels: Same spectrum of 2MASS J1534—2 9 52AB plotted as a thick black 
line. Spectra of T dwarf spectral standards from iBurgasser et al.l (|2006bl ) are plotted as colored 
lines. The spectra have been normalized by their peak flux. 
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Fig. 3. — Ne ar-IR colors of 2MA SS J1534— 2952AB co mpared with nearby single late-L and 
T dwarfs fron i iKnapp et al.l (12004 ) and Chiu et al.l (120061) a.nd ind i vidual c omponents of resolve d 
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Liu et al.l (|2006bl ) . The photometry errors are comparable to or smaller than the size of the plotting 
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being rounded down (e.g., T3.5 is labeled as "3") , and objects of the same subclass plotted in the 
same color. The late-L dwarfs (classified on the iGeballe et al.l l2002l scheme) are plotted as bare 
numbers. The T dwarfs (on the IBurgasser et al.ll2006bl scheme) are plotted as circumscribed num- 



bers, with squares for integer subclasses (e.g., T3) and diamonds for half subclasses (e.g., T3.5). 
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Fig. 4. — {CH4S — H) color versus near-IR spectral type for T dwarfs. The plotting symbols are 
the synthesized colors from published low-resolution near-IR spectra of T dwarfs, excluding objects 
that are spectrally peculiar and known binaries. (See § 13.11 for references.) The lines represent 
2nd-order polynomial fits to the data, with the solid line being the spectral type as a function of 
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Fig. 5. — Resulting MCMC probability distributions of orbital parameters: semi-major axis (a), 
orbital period (P), eccentricity (e), inclination (i), epoch of periastron (Tq), PA of the ascending 
node (0), and argument of periastron (w). Each histogram is shaded to indicate the 68.3% and 
95.5% confidence regions, which correspond to 1 and 2a for a normal distribution, and the solid 
vertical line represents the median value. 
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Fig. 6. — Results from MCMC determination of the orbital period and eccentricity for 
2MASS J1534— 2952AB, illustrating the degeneracy between the two parameters. The central plot 
shows all the values in the MCMC chain. Two branches of possible orbits are seen, a short-period 
{P < 20 yr) branch and a long-period (P > 20 yr) one. About 98% of the MCMC chain steps 
are in the short-period branch. The top and side plots show the resulting probability distributions 
of P and e. Each histogram is shaded to indicate the 68.3% and 95.5% confidence limits, which 
correspond to la and 2a for a normal distribution, and the dashed vertical lines represent the 
median values. 
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Fig. 7. — Results from MCMC determination of orbital period and the argument of periastron 
illustrating the degeneracy between these two parameters. The locus illustrates the two general 
branches of possible orbits, with short-period orbits having just completed apoastron and long- 
period orbits having just completed periastron. See Figure E] for further explanation. 
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further explanation. The dashed lines represent the ±3(7 confidence intervals of the mass probability 
distribution. 
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Fig. 9.— Keck LGS (red dot) and HST (blue square) measurements of 2MASS J1534-2952AB 
along with three representative orbits from the MCMC fitting. Measurement errors are comparable 
to or smaller than the colored plotting symbols. The solid line shows the best-fitting orbit with a 
15.2-year period, and the two dashed orbits represent alternative long-period (20 year) and short- 
period (12 year) solutions, chosen to represent the plausible range in period. The three orbits have 
reduced chi-square values of 0.9, 1.0, and 1.1. The empty circles are the location of the secondary in 
future years as predicted by the 15.2-year period orbit. The image insets are 1" on a side, displayed 
with a square-root stretch. (The HST/WFPC2 image cutout has been rotated so that North is up. 
See Figure [U for the most accurate representation of the original data.) 



49 



250 







■ HST 




• <eck ^y^^'^'^'^'*-^. 


^ 200 




w 


^X'' 'I 


CO 


//•■ \ 


E 


r' V 


^ 150 


'A \ 


c 


//■■ \ 


o 


a/^ \ 


-I— ■ 


/ \ 


2 100 


'/ \ 


CO 


£ \ 


Q. 


/ \ / 





I \ /A 


^^ 50 


'■.f \f 







w 




1^ 

V -2 


"V '^■*- 1- 4- ^f ■■■■'' '<■ '" 


t. ^^..:;;r..i-.-.-«,,,-^«^ ^.^ , 


O 




20 


00 2002 2004 2006 2008 2010 




Year 




0.0 
O -0.5 
-1.0 



2000 2002 2004 2006 2008 2010 
Year 



Fig. 10.— Keck LGS and RST measurements of 2MASS J1534-2952AB's separation {left) and PA 
(right), along with the three orbits plotted in Figure [9j The measurement errors are comparable 
to or smaller than the plotting symbols, except for the three data points shown in the PA plot. 
The bottom panel of each plot shows the difference of the observed astrometry and the best-fitting 
(15.2-year period) orbit. 
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Fig. 11. — Total mass distribution from our MCMC analysis. The histogram is shaded to indicate 
the derived 68.3%, 95.5%, and 99.7% confidence regions, which correspond to la, 2a, and 3a for a 
normal distribution. The solid vertical line represents the median value of 0.0556 Mq. The standard 
deviation of the distribution is 0.0018 Mq. The wider, unshaded histogram shows the final mass 
distribution, after accounting for the additional 4.9% error due primarily to the uncertainty in the 
binary's parallax — the result is essentially gaussian with a standard deviation of 0.003 Mq. The 
confidence limits for both distributions are given in Table [H 
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Fig. 12. — Total mass distribution from the MCMC analysis using four different priors: a distribu- 
tion that is flat in the logarithm of the period and se mi-major axis (our de f ault a ssumption); the 
log- normal period distribution for solar- type stars from lDuaueiinov &: Mavorl (Il99ll ): the log- normal 
semi-major distribution for field u ltracool binaries fromlAllenI ()2007l ): and the truncated log- normal 



semi-major axis distribution from lMaxted fc JeffriesI (J2005l ). The median value is indicated by the 
filled circle and the two sets of errors bars indicate the 68.3% and 95.5% confidence regions, which 
correspond to la and 2a for a normal distribution. The dashed curves at the bottom sho w the 
result ing mass distributions after accounting for the uncertainty in the parallax. While the lAllen 
(J2007l ) prior favors slightly higher masses, all four priors give very consistent results. 
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Fig. 13. — Determination of the age of the 2MASS J1534— 2952AB system from evolutionary models 
using the ob served magnitudes /lu minosity and total mass. T wo sets of models ar e employed: the 
models from lBurrows et al.l ( 19971 ) and the COND models from lBaraffe et al.l (120031 ). The horizontal 
grey bars indicate our measured 1, 2, and 3a constraints on the total mass, and the vertical lines 
show the corresponding ages derived from the Burrows et al. models. Both sets of models produce 
consistent results, given in Table [5j The COND models provide predictions for both the absolute 
magnitudes and bolometric luminosity, so all of these are shown. (The two sets of models predict 
very similar L^oi results, so the red and black model lines are indistinguishable on this plot.) Note 
that these plotted curves are computed for an object with the fluxes of 2MASS J1534— 2952AB and 
are not generally applicable to other binaries. See § 14.21 for details. 
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Fig. 14. — Color-magnitude diagrams for the individual components of 2MASS J1534— 2952AB 
compared to the Lyon COND models ([Baraffe et al.ll2003l ). on the CIT photometric system. The 
models corresponding to the individual component masses and their ztlcr range are plotted as solid 
lines and dotted lines, respectively. The numbered boxes indicate the logarithm of the model age 
in Gyr. The error bars on the absolute J-band magnitude is smaller than the plotting symbol. The 
models do not match the data very well, which can be ascribed to the incomplete opacities in the 
model atmospheres. 
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Fig. 15. — H-R diagram showing the individual components of 2MASS J1534— 2952AB compared 
to evolutionary tracks from the Tucson and Lyon groups. The plotted symbols represent the data, 
with the effective temperatures determined from the spectral types and model atmosphere studies 
of field T dwarfs (Equation [7]) . The uncertainity in the Lboi measurements are comparable to the 
vertical extent of the plotting symbols. The solid line represents the median mass values determined 
for the individual components, and the dotted lines show the itlo" ranges. The embedded small 
squares demarcate logarithms of a set of ages; the ages for the Burrows models are labelled, and 
the (unlabeled) ages for the Lyon models have very similar positions along the plotted line. There 
is a modest systematic disagreement between the evolutionary tracks and observations. 
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Fig. 16. — Contour plots showing the uncertainties in determining T^g and log{g) using mass 
benchmarks (brown dw arfs with dynamical in asses) and age benchmarks (brown dwarfs with known 
ages), derived from the I Bur rows et al.l (|200ll ) analytic approximations to substellar evolution. The 
X- and y-axes give the fractional error in the measurements of mass, age, and -Lboi, and the contours 
show the resulting fractional uncertainties in T^n and log{g) spaced from 1% to 40%. 



Table 1. Keck LGS AO Observations 



Date 
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Airmass 


FWHM 


Strehl ratio 


Separation'^ 


Position angle'' 


Amag 


(UT) 






(mas) 




(mas) 




(dog) 




2005-May-Ol 


J 


1.66 


102 ± 12 


0.020 ± 0.002 


211.3 ± 1.5 


(1.5) 


14.1 ± 0.3 (0.3) 


0.163 ±0.014 




H 


1.63 


86 ±6 


0.047 ± 0.006 


211.7 ± 0.8 


(0.8) 


13.86 ± 0.15 (0.13) 


0.286 ± 0.011 




K' 


1.61 


88 ±6 


0.101 ± 0.012 


212.4 ± 1.1 


(1.0) 


14.0 ± 0.2 (0.2) 


0.278 ± 0.021 


2006-May-05 


Ks 


1.56 


64 ±3 


0.210 ± 0.014 


190.6 ± 0.3 


(0.2) 


15.43 ± 0.12 (0.09) 


0.282 ± 0.010 


2007-Mar-26 


K 


1.56 


82 ±3 


0.151 ± 0.016 


158.0 ± 0.6 


(0.6) 


17.5 ± 0.2 (0.19) 


0.287 ± 0.012 


2007-Apr-22 


Ks 


1.57 


67 ±5 


0.20 ± 0.03 


153.7 ± 0.4 


(0.3) 


17.53 ±0.13 (0.10) 


0.269 ± 0.010 


2008-Jan-15 


Ks 


2.05 


100 ± 3 


0.074 ± 0.002 


114.4 ± 1.1 


(1.1) 


21.5 ± 0.9 (0.9) 


0.27 ± 0.06 


2008-Apr-Ol 


Ks 


1.55 


87 ± 4 


0.095 ± 0.018 


102.5 ± 0.7 


(0.7) 


21.1 ± 0.7 (0.7) 


0.25 ± 0.04 




CHiS 


1.58 


78 ± 7 


0.048 ± 0.018 


102.0 ± 0.4 


(0.4) 


20.4 ± 1.5 (1.5) 


0.21 ± 0.04 



^All photometry on the MKO system. 

The tabulated errors are computed by appropriately combining in quadrature: (1) the instrumental measurements from 
fitting the images of the binary and (2) the overall uncertainties in the NIRC2 pixel scale and orientation. The errors in 
parentheses represent the instrumental errors alone. See § 2 for details. 
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Table 2. HST Observations 



Date Instrument Filter Separation^ Position angle*^ 

(UT) (mas) (mas) 



2000-Aug-18 
2006-Jan-19 
2006- Apr- 11 



WFPC2 
ACS 
ACS 



F8UW 
F8UW 
F8UW 



62.8 ± 1.2 
199.0 ± 1.1 
191.2 ± 1.1 



357.1 ± 0.8 
14.5 ± 0.6 
15.5 ± 0.4 



Amae 



0.30 ± 0.05 
0.28 ± 0.06 
0.30 ± 0.04 



'^The tabulated errors are dominated by the uncertainties in fitting the binary im- 
ages, which are much larger than the errors in the overall astrometric calibration of 
WFPC2 and ACS. See S 2 for details. 
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Table 3. Resolved Properties of 2MASS J1534-2952AB'' 

Property 2MASS J1534-2952A 2MASS J1534-2952B 

F8UW - J (mags) 4.95 ± 0.04 5.10± 0.04 

J- H (mags) -0.08 ± 0.04 -0.21 ± 0.04 

CH4S - H (mags)= -0.30 ± 0.05 -0.37 ± 0.05 

H - K {TCiSigs) -0.17 ± 0.04 -0.17 ± 0.04 

J -K (mags) -0.25 ± 0.04 -0.38 ± 0.04 

Mfsiaw (mags) 19.57 ± 0.04 19.87 ± 0.05 

Mj (mags) 14.61 ± 0.05 14.77 ± 0.05 

M// (mags) 14.69 ± 0.05 14.98 ± 0.05 

Mk (mags) 14.86 ± 0.05 15.15 ± 0.05 

Estimated spectral type"" T5.0 ± 0.5 T5.5 ± 0.5 

log(Lboi/^0)'' -5.015 ±0.019 -5.093 ±0.019 

'^All infrared photometry on the MKO photometric system. 

'^Based on the iBurgasser et al.l ( 2006bl ) near-IR classification scheme. 



■^The difference in the {CH4S — H) color of the two components is 0.07 ± 
0.02 mags, i.e., better constrained than the quadrature sum of the measure- 
ment errors tabulated here. See ^ 13.11 Similarly, the difference in log(Lboi/^0) 
is 0.078 ± 0.016 dex, since this quantity is independent of the distance uncer- 
tainty. 
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Table 4. Derived Orbital Parameters for 2MASS J1534-2952AB 







MCMC 




ORBIT 




Median 


68.3% c.l. 


95.5% c.l. 




Time of periastron Tq - 2400000.5 (MJD) 


55960^ 


-240, 210 


-740, 450 


56024 ± 347 


Orbital period P (yr) 


15.1 


-1.6, 2.3 


-3.1, 5.1 


15.2 ± 2.6 


Semi-major axis a (mas) 


171 


-13, 19 


-27, 41 


172 ± 22 


Semi-major axis a (AU)*^ 


2.3 


-0.2, 0.3 


-0.4, 0.6 


2.3 ± 0.3 


Inclination i (°) 


84.3 


-0.6, 0.6 


-1.7, 1.0 


84.3 ± 0.8 


Eccentricity e 


0.25 


-0.13, 0.11 


-0.20, 0.25 


0.24 ± 0.16 


PA of the ascending node 17 (°) 


13.0 


-0.3, 0.3 


-0.9, 0.5 


13.0 ± 0.4 


Argument of periastron uj (°) 


179 


-14, 6 


-83, 11 


178 ± 10 


Total mass (Mq): fitted 


0.0556 


-0.0017, 0.0019 


-0.004, 0.004 


0.056 ± 0.004 


Total mass {Mq): final 


0.056 


-0.003, 0.003 


-0.006, 0.007 


0.056 ± 0.005 



^February 3, 2012 UT. 

"^Includes the uncertainty in the parallax and pixel scale. 



Note. — Median values and confidence limits for orbital parameters derived from our default MCMC 
fitting, which uses a prior distribution flat in log(P) and log(a). The "fitted total mass" represents the 
direct MCMC results from fitting the observed orbital motion of the two components. The "final total 
mass" includes the additional 4.9% error from the uncertainties in the parallax and the Kcck/NIRC2 pixel 
scale; the former is R:!l5x larger than the latter. The final m ass distribution is esse ntially gaussian. The 
rightmost column gives the results from the ORBIT routine bv lForveille et al.l (J1999I ). (See ? 13.3.21 ) 



Table 5. Evolutionary Model-Derived Properties of 2MASS J1534-2952AB 



Property 



Tucson models 
Component A Component B 



Lyon models 
Component A Component B 



"Average" 



Component A 



Component B 



log(age) 

Radius {Rq) 0.0997 

Teff (K) 1019 

log(5) (cgs) 4.90 

ATeff (K) 

Mass ratio Mb /Ma 

Mass (Mq) 0.0287 



8.86 

+0.0012(0.0025) 
-0.0012(0.0024) 

+ 16(32) 
-15(31) 

+0.04(0.07) 
-0.04(0.07) 



+0.04(0.09) 
-0.04(0.09) 

0.1004 



970 



-0.0012(0.0024) 
-0.0012(0.0024) 

-15(31) 
-15(29) 



50 
0.937 

+0.0017(0.0036) 
-0.0016(0.0032) 



A oc +0.04(0.07) 
^■°" -0.04(0.07) 

+6(10) 
-10(19) 

+0.012(0.024) 
-0.007(0.012) 

+0.0016(0.0033) 
0.0015(0.0030) 



0.0270 



;.92 



+0.04(0.09) 
-0.04(0.09) 



n 0078 +0'0010(0.0021) „ „„„ . +0.0011(0.0022) 

u.uy(o -0.0010(0.0022) u.uyoy; _o. 0012(0.0025) 



1034 +i'^(32) 



Qor +15(32) 
^°^ -15(29) 



4.92 



HO. 04(0. 07) 
-0.04(0.07) 



+0.04(0.08) 
-0.04(0.07) 



50 



+6(10) 
-10(19) 



0.934 



+0.012(0.025) 
-0.008(0.013) 



n 0988 +0.0017(0.0035) p. f,r)f.Q +0.0017(0.0034) 

U.UZOO -0.0016(0.0031) U.U^D» _o. 0016(0. 0030) 



8.89 ±0.05(0.10) 
0.0986 ± 0.0015(0.0028) 0.0993 ± 0.0017(0.0031) 
1028 ± 17(35) 978 ± 17(34) 

4.91 ± 0.04(0.07) 4.87 ±0.04(0.07) 



50 



+6(10) 
-10(19) 



0.936 



+0.012(0.024) 



-0.008(0.013) 

0.0287 ± 0.0016(0.0033) 0.0269 ± 0.0016(0.0032) 



I 

O 

Note. — Median values of physical parameters derived for the two components of 2MASS J1534— 2952AB from the evolutionary models along with the (j)8(95%) confidence 
limits, as described in § 14.21 By construction, the ages of the two components are identical for a given set of models, since the analysis assumes the system is coeval, and thus a 
single model-derived age is listed for both components. Note that the two sets of models predict nearly identical mass ratios so the two sets of individual masses are the same. 



